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Preambulo 



La presente memoria esta dedicada al estudio de los efectos inducidos por la 
masa de los quarks en observables relacionados con la desintegracion a tres jets 
del boson de gauge Z y su aplicacion experimental en LEP. 

Para la mayor parte de los observables estudiados en LEP el efecto de la 
masa de los quarks puede despreciarse puesto que esta siempre aparece como el 
cociente m 2 /m z , donde m q es la masa del quark y es la masa del boson 
Z, lo cual supone correcciones muy pequenas incluso para el quark mas pesado 
producido en LEP, el quark 6, muy por debajo de la precision experimental 
actualmente accesible. 

No obstante, aunque este ultimo argumento sea cierto para secciones eficaces 
totales no ocurre asi cuando estudiamos observables que, aparte de la masa del 
quark, dependan de variables adicionales como por ejemplo secciones eficaces a 
n-jets. En tal caso, puesto que ponemos en juego una nueva escala de energfas 
E c = mz^/Vc, donde y c es el parametro que define la multiplicidad del jet, 
aparecen contribuciones del tipo m 2 /E 2 = (m q /m z ) 2 jy c que, para valores de 
y c lo suficientemente pequenos, podrfan incrementar considerablcmente el efecto 
de la masa de los quarks y permitir su estudio en LEP. Asf mismo, el efecto 
de la masa de los quarks podrfa verse favorecido por logaritmos de la masa, 
m q \og(m q /m z ), provenientes de la integration sobre espacio fasico. 

Dos son las razones que motivan nuestro estudio. En primer lugar, el error 
sistematico mas importante en la medida de la constante de acoplamiento fuerte 
a b s (mz) a partir del cociente entre las anchuras de desintegracion del boson Z a 



dos y tres jets en la production de bb en LEP [ 130 | precede de la incertidumbre 
debida al desconocimiento de los efectos de la masa del quark b. Una mejor 
comprension de estos efectos contribuiria obstensiblemente a una mejor medida 
de la constante de acoplamiento fuerte a h s (mz)- 

En segundo lugar, asumiendo universalidad de sabor para las interacciones 
fuertes, dicho estudio podn'a permitir por primera vez una medida experimental 
de la masa del quark b a partir de los datos de LEP. Debido a que los quarks no 
aparecen en la naturaleza como partfculas libres el estudio de su masa presenta 
serias dificultades teoricas. De hecho, la masa de los quarks debe considerarse 
como una constante de acoplamiento mas. La masa de los quarks pesados, 
como el b y el c, puede extraerse a bajas energfas a partir del espectro del 
bottomium y del charmonium con Reglas de Suma de QCD y calculos en el 
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reti'culo. Una medida de la masa del quark b a altas energi'as, como por ejemplo 
en LEP, presentaria la ventaja de permitir una determination de la masa del 
quark b a una escala de energi'as muy por encima de su umbral de production, a 
diferencia de lo que ocurre en los dos metodos anteriormente descritos. Es mas, 
dicha medida permitirfa por primera vez comprobar como la masa del quark b 
evoluciona segun predice el Grupo de Renormalizacion desde escalas del order 
de la masa del quark misma, fi = rrib, hasta altas energfas, [i = mz, del mismo 
modo que fue posible hacerlo para la constante de acoplamiento fuerte a s (p) 
y constituirfa una nueva confirmation de QCD como teoria para describir las 
inter acciones fuertes. 

En el Capi'tulo 1 comenzamos revisando las distintas definitiones de masa, 
repasamos cuales son las Ecuaciones del Grupo de Renormalizacion en QCD, 
analizamos como conectar los parametros de una teoria con los de su teoria 
efectiva a bajas energi'as mediante las adecuadas Condiciones de Conexion y re- 
alizamos una pequena recopilacion de las determinaciones mas recientes que de 
las masas de todos los quarks han sido realizadas a partir de Rcglas dc Suma de 
QCD, calculos en el reti'culo y Teon'a de Perturbaciones Quirales. Finalmente, 
evolucionamos todas las masas mediante las Ecuaciones del Grupo de Renor- 
malizacion hasta la escala de energfas de la masa del boson Z poniendo especial 
incapie en la masa del quark b cuyos efectos pretendemos estudiar en profundidad 
en el resto de esta tesis. 

En el Capi'tulo 2 justificamos por que es posible medir la masa del quark b 
en LEP a partir de observables a tres jets, comparamos el comportamiento de la 
anchura de desintegracion inclusiva del boson Z expresada en funcion de las dis- 
tintas definitiones de masa, definimos cuales son los algoritmos de reconstruction 
de jets, en particular los cuatro sobre los cuales basaremos nuestro analisis: EM, 
JADE, E y DURHAM. Finalmente, analizamos a primer order en la constante de 
acoplamiento fuerte algunos observables a tres jets y proporcionamos funciones 
sencillas que parametrizan su comportamiento en funcion de la masa del quark 
y del parametro que define la multiplicidad del jet. 

Puesto que el primer order no permite distinguir entre las posibles defini- 
tiones de masa y puesto que para el quark b la diferencia cntre cllas es numcricamente 
muy importante, en los Capftulos 3, 4 y 5 nos centramos en el analisis del orden 
siguiente. En el Capi'tulo 3 presentamos y clasificamos las distintas amplitudes 
de transition que debemos calcular a este orden. La principal dificultad de dicho 
calculo radica en la aparicion de divergencias infrarrojas debido a la presen- 
cia de parti'culas sin masa como son los gluones. En el Capftulo 4 analizamos 
el comportamiento infrarrojo de estas amplitudes de transition, la integramos 
analfticamente en la region de espacio fasico que conticne las divergencias y fi- 
nalmente mostramos como dichas divergencias se cancelan cuando sumamos las 
contribuciones procedentes de las correcciones virtuales y reales. En el Capi'tulo 
5 presentamos los resultados de la integration numerica de las partes finitas asf 
como ajustes sencillos de estos resultados para facilitar su manejo, exploramos a 
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segundo orden los observables a tres jets que habfamos estudiado en el Capftulo 
2 y discutimos su comportamiento en funcion de la escala. Finalmente, en los 
apendices, recopilamos algunas de las funciones necesarias para el calculo que 
hemos realizado: integrales a un loop, reduction de Passarino-Veltman, el prob- 
lema de 75 en Regularizacion Dimensional, espacio fasico en _D-dimensiones, etc. 

El resto de esta tesis esta escrita en ingles para cumplir con la normativa 
vigente sobre el Doctorado Europeo. 
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Preamble 



This thesis is devoted to the study of the effects induced by the quark mass 
in some three-jet observables related to the decay of the Z gauge boson and its 
experimental application in LEP. 

Quark masses can be neglected for many observables at LEP because usually 
they appear as the ratio m^/ml, where m q is the quark mass and mz is the 
Z-boson mass. Even for the heaviest quark produced at LEP, the &-quark, these 
corrections are very small and remain bellow the LEP experimental precision. 

While this argument is correct for total cross sections, it is not completely 
true for quantities that depend on other variables. In particular n-jet cross 
sections. In this case we introduce a new scale E c — raz^/Wc, where y c is the 
jet-resolution parameter that defines the jet multiplicity, and for small values of 
y c there could be contributions like m^/E 2 = (m q /m z ) 2 /y c which could enhance 
the quark mass effect considerably allowing its study at LEP. Furthermore, the 
quark mass effect could be favoured by logarithms of the mass, m q log (m^/m|), 
coming from phase space integration. 

Our motivation is twofold. First, it has been shown |130| that the biggest 
systematic error in the measurement of a b s (mz) (a s obtained from 66-production 
at LEP from the ratio of three to two jets) comes from the uncertainties in the 
estimate of the quark mass effects. A better understanding of such effects will 
contribute to a better determination of the strong coupling constant a b s (mz)- 

Second, by assuming flavour universality of the strong interaction we study 
the possibility of measuring the bottom quark mass, to?,, from LEP data. A 
precise theoretical framework is needed for the study of the quark mass effects in 
physical observables because quarks are not free particles. In fact, quark masses 
have to be treated more like coupling constants. The heavy quark masses, like 
the b and the c-quark masses, can be extracted at low energies from the botto- 
mia and the charmonia spectrum from QCD Sum Rules and Lattice calculations. 
Nevertheless, a possible measurement of the bottom quark mass at high energies, 
like in LEP, would present the advantage of determining the bottom quark mass 
far from threshold in contrast to the other methods described above. Further- 
more, such a measurement would allow to test for the first time how the bottom 
quark mass evolve following the Renormalization Group prediction from scales 
of the order of the quark mass itself, y, — mb, to high energy scales, = mz, 
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in the same way it was possible with the strong coupling constant and would 
provide a new test of QCD as a good theory describing the strong interaction. 

In Chapter 1 we start by reviewing the different theoretical mass definitions, 
we solve the QCD Renormalization Group Equations, we analyze how to connect 
the parameters of a theory with the parameters of its low energy effective theory 
through the adequate Matching Conditions, in particular how to pass a heavy 
quark threshold and we review the most recent determinations of all the quark 
masses form QCD Sum Rules, Lattice, and Chiral Perturbation Theory. Finally 
we run all the quark masses to the Z-boson mass scale and we focus our study 
in the bottom quark mass which we will study in the rest of this thesis. 

In Chapter 2 we justify the possibility of extracting the bottom quark mass in 
LEP from some three-jet observables, we compare the behaviour of the Z-boson 
inclusive decay width expressed in terms of the different quark mass definitions, 
we define what is a jet clustering algorithm, in particular the four on which we 
will work: EM, JADE, E and DURHAM. Finally, we analyze at first order in the 
strong coupling constant some three-jet observables and we give simple functions 
parametrizing their behaviour in terms of the quark mass and the jet resolution 
parameter. 

Since the first order does not allow to distinguish which mass definition 
should we use in the theoretical expressions and since for the bottom quark 
mass the difference among the possible mass definitions is quite important, in 
Chapters 3, 4 and 5 we focus in the analysis at the next order. In Chapter 3 
we present and classify the transition amplitudes we should compute at second 
order in the strong coupling constant. The main difficulty of such calculation 
is the appearance, in addition to renormalized UV divergences, of infrared (IR) 
singularities since we are dealing with massless particles like gluons. In chapter 
4 we analyze the infrared behaviour of our transition amplitudes, we integrate 
them analytically in the phase space region containing the infrared singularities 
and finally we show how the infrared divergences cancel when we sum up the 
contributions coming from the real and the one-loop virtual corrections. In 
Chapter 5 we present the results of the numerical integration over the finite 
parts and we give simple fits to these results in order to facilitate their handling, 
we explore at second order the three-jet observables we have studied in Chapter 
2 and we discuss their scale-dependent behaviour. Finally, in the appendices, 
we collect some of the functions needed for the calculation we have performed: 
one-loop integrals, Passarino-Vcltman reduction, the 75 problem in Dimensional 
Regularization, phase space in D-dimensions, etc. 

This thesis fulfils the European Ph.D. conditions. 



Contents 



Preambulo 



Preamble 



Chapter 1. The quark masses 



Chapt 


er 2. Three-jet observables at leading order 




1. 


The inclusive decay rate Z 


-> bb 






2. 


Jet clustering algorithms 










3. 


Two- and three-jet event rates 






4. 


Discussion and conclusions 







Chapter 3. Transition amplitudes at next-to- leading order 



1. Virtual corrections 

2. Emission of two real gluons 

3. Emission of four quarks 



Chapter 4. Infrared cancellations 



1. Soft divergences 



2. Collinear divergences 



Chapter 5. 


Numerical results and conclusions 






1. Tests 










2. Massive three-jet decay rate at next-to- leading order 




3. Discussion and conclusions 





Appendix A. Loop integrals 

One- and two-point functions] 



2. Three-point functions 

3. Four-point functions 



4. Explicit calculation of C03 

5. Passarino-Vcltman reduction 



6. The problem of 75 in Z?-dimcnsions 



Appendix B. Phase space in D = 4 — 2e dimensions 




1. System 1-3 





11 

13 
15 
18 
25 

27 
28 
30 
33 

39 
39 
44 

■51 
53 
57 
GO 

67 
67 
G8 
71 
72 
72 
74 

75 
75 



viii 



2. System 3-4 



CONTENTS 



Bibliography 



CHAPTER 1 



The quark masses 



GUT and SUSY theories predict some relations among the fermion masses 
(or more properly among Yukawa couplings) at the unification scale, Mqut ~ 
10 16 (GeU). For instance, in SU(5) we have the usual lepton-bottom quarks 
unification, hi, — h T , h s — h^, hd = h e , or the modified Gcorgi-Jarlskog relation, 
hb = h T , h s = h/j,/3, hd = 3/i e , while for 5*0(10) typically we get unification of 
the third family, ht = hb = h T = h Vr . This together with the RGE's provide a 
powerful tool for predicting quark masses at low energies. 

This chapter is not devoted to Unification. Rather, we try to make a short 
review of the most recent determinations of the quark masses and to calculate 
the running until mz, the mass of the Z-boson. Why mz?. For model building 
purposes it is a good idea to have a reference scale and extensions of the Standard 
Model appear above mz- Furthermore, the strong coupling constant a s is really 
strong below mz and then special care has to be taken on the matching in passing 
a heavy quark threshold and on the running. 

To define what is the mass of a quark is not an easy task because quarks 
are not free particles. For leptons it is clear that the physical mass is the pole 
of the propagator. Quark masses, however, have to be treated more like cou- 
pling constants. We can find in the literature several quark mass definitions: the 
Euclidean mass, Me(p 2 — — M 2 ), defined as the mass renormalized at the Eu- 
clidean point p 2 = —M 2 . It is gauge dependent but softly dependent on Aqcd- 
Nevertheless, it is not used anymore in the most recent works. The perturbative 
pole mass [126| and the running mass are the two most commonly used quark 
mass definitions. The "perturbative" pole mass, M(p 2 = M 2 ), is defined as the 
pole of the renormalized quark propagator in a strictly perturbative sense. It is 
gauge invariant and scheme independent. However, it suffers from renormalon 
ambiguities. The running mass, m(fi), the renormalized mass in the MS scheme 
or its corresponding Yukawa coupling related to it through the vev of the Higgs, 
m(/i) = v(/j,)h(lJ,), is well defined and since it is a true short distance parameter 
it has become the preferred mass definition in the last years. 

For the light quarks, up, down and strange, chiral perturbation theory 51, 



52 



, 61 -63 1 provides a powerful tool for determining rcnormalization group in- 
variant quark mass ratios. The absolute values, usually the running mass at 
l(GeU), can be extracted from different QCD Sum Rules §, g, |o|, |l0^| 

or for the strange quark mass from lattice Q]. For the heavy quarks, bottom 



l 



2 



1. THE QUARK MASSES 



Table 1 . Recent determinations of the light quark masses from 
second order xPT, QCD Sum Rules and lattice. 



Gasser 
Leutwyler 


XPT 


0( P 4 ) 


m d - m M 2m a = 2M x 1Q -3 
m s — m m s + m 

m s /m = 25.7 ± 2.6 


Donoghue 
et al. 


XPT 


0( P A ) 


m d -m u 2m _ n n x 1Q -3 
m s — m m s + rh 

m s /m = 31. 



Bijnens 
Prades 
de Rafael 


FESR 
Laplace SR 
(pseudo) 


NNLO 


{fh u +m d )(lGeV) 
= 12.0 ±2.5 


a s {m z ) = 0.117(5) 


loffe 
et al. 


Isospin viol, 
in QCD SR 


NNLO 


{fhd - fh u )(lGeV) 
= 3.0 ± 1.0 


A = 150. 



Narison 


r-like SR 


NNLO 
NLO 


m s (lGeV) 
fh s (lGeV) 


= 197(29) 
= 222(22) 


a s {m z ) = 0.118(6) 


Jamin 


QSSR 


NNLO 


fh s {lGeV) 


= 189(32) 


a s (m z ) = 0.118(6) 


Miinz 


(scalar) 










Chetyrkin 


QCD SR 


NNLO 


fh 8 {lGeV) 


= 171(15) 


a s (m z ) = 0.117(5) 


et al. 












Allton 


quenched 


NLO 


fh s (2GeV) 


= 128(18) 


A 5 = 240. ± 90. 


et al. 


Lattice 











and charm, we can deal either with QCD Sum Rules [48-50, 
or lattice calculations [42-44,53, 



68 



109 111 



127] 



The bulk of this thesis is devoted to 
explore the possibility of extracting the bottom quark mass from jet physics at 
LEP [^5, 59, SO, 12C[ |. For the top quark we have the recent measurements from 
CDF and D0 at FERMILAB ^ |65|, |98|, pi] that we will identify with the 
pole mass. 

We have summarized in tables |l] and ^| all of these recent quark mass determi- 
nations. Of course the final result depends on the strong gauge coupling constant 
used in the analysis, for this reason we quote it too. In the running we will take 



the world average [22| strong coupling constant value a s (mz) = 0.118 ± 0.006 
for masses obtained from QCD Sum Rules but for lattice masses we will run with 
the lattice |44] result a[°'(mz) 



0.115±0.002. These values are consistent with 
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Table 2. Recent determinations of the heavy quark masses 
from QCD Sum Rules, lattice and FERMILAB. 



Narison 


QSSR 
T 


NLO 


m b {M b ) = 4.23(4) 
m c (M c ) = 1.23+g 
M b = 4.62(2) 
M c = 1.42(3) 


a s {m z ) = 0.118(6) 


Narison 


non-rel 
Laplc.SR 


NLO 


M™ = 4.69+gj 
M? R = 1.45+gj 


a s (m z ) = 0.118(6) 


Dominguez 
et al. 


rel, non-rel 
Laplc.SR 
J/*, T 


LO 

l/m 2 q 


M b = 4.70(7) 
M c = 1.46(7) 


A 4 = 200 - 300 
A 5 = 100 - 200 


Titard 
Yndurain 


potential 




m b (m b ) = 4.397+I33) 
m c (m c ) = 1.306+ ( ( 35) 


a s {m z ) = 0.117(5) 


Neubert 


QCD SR 


NLO 

l/m q 


M 6 = 4.71(7) 
M c = 1.30(12) 





Crisafulli 
et al. 


Lattice 
in B-meson 


m b (fh b ) = 4.17(6) 




Gimenez 
et al. 


Lattice 
in B-meson 


m b (fn b ) = 4.15(20) 




Davies 
et al. 


NRQCD + leading 
rel and Lattice 
spacing, bb 


M b = 5.0(2) 
fh b (M b ) = 4.0(1) 


0^ = 0.115(2) 


El-Khadra 
Mertens 


Fermilab action 
in quenched Lat 


M c = 1.5(2) 





CDF 


M t = 176.8 ± AA{stat) ± 4.8(sys) 


mean 


D0 


M t = 169. ± 8.(stat) ± 8.(sys) 


M t = 175. ±6. 



almost all the references. For those that differ an updat e is n eeded but this is be- 
yond the goals of this review. For instance, S. Narison [109] makes two different 



determinations for the bottom and the charm quark masses. In the first, and for 
the first time, he gets directly the running mass avoiding then the renormalon 
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ambiguity associated with the pole mass. The second one, from non-relativistic 



Laplace Sum Rules, is in fact an update of the work of Dominguez et al. [48-50 

The 0(a 2 ) strong correction to the relation between the perturbative pole 
mass, M, and the running mass, fh(fJ,), was calculated in U~ 



m(M) 3 

where K t ~ 10.95 for the top quark, Kj, ~ 12.4 for the bottom and K c ~ 13.3 for 
the charm. As pointed out by S. Narison [109] cq. (fhl]) is consistent with three 



loops running but for two loops running we can drop the 0(a 2 s ) term. Recently, 
the electroweak correction to the relation between the perturbative pole mass 



and the Yukawa coupling has been calculated [28, 75 1. However, this correction 
is small, for instance for the top quark it is less than 0.5% in the SM for a mass 
of the Higgs lower than 600(GeV) and at most 3.% for M H ~ l(TeV), and for 
consistency one has to include it only if two loop electroweak running is done. 

Instead of expressing the solution of the QCD renormalization group equa- 
tions for the strong gauge coupling constant and the quark masses in terms of 
Aqcd we perform an expansion in the strong coupling constant at one loop 



12l| . At three loops we get 

a s ( M ) = a«( M )(l + Cl ( M )aW(/x) + c 2 ( M )(aW( M )) 2 ), (1.2) 
(/i) = m( 1 )( M )(l + d 1 ( M )ai 1 )( Ai )+d 2 ( A1 )(aW( A1 )) 2 ), (1.3) 



nil 



where a^ 1 (/i) and mP^(ji) are the one loop solutions 

= i. "^/,, ' ™ (1) 0*) = m(Mo)^(M)- 2 ^ //3 °, (1.4) 
1 + a s (fi )Pot 

with t = l/(47r) log/x 2 //itg, K{fJt) the ratio K(fj,) = a s {ii$) / cSp^), and 



c 2 ^)=bl\ogK{Li) [logK(fi) - 1] 
270 



di(n) 



A) 



i(bi 



{b\ b 2 ) [1 - KM] 
gi)[l-K(ji)]+b 1 log K (/*)], 



(1.5) 



7o 



W = f\lPo(b 2 - K) + 270(^1 - 5i)1 [1 - K(ji)]' 
+ Po(g2-b igi )[l-K 2 (n)] 

+ 4 70 6i(6i - 31) [1 - K{ji)\ - 2/3 6i3i + 6?(/3o + 2 7 o) logi^) 



log K (n) , 
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where 

u Pi , fo 7i 72 n R , 

bi = —5-, 02 = ,. , 2n , 9i = , 92 = -5—, (1.6) 

4tt/3 (4tt) 2 /3 tt7o n 2 j 

are the ratios of the well known beta and gamma functions in the MS scheme 

00 = 11 -^V>, 70 = 2, 

38 1 / 5033 325 - 

01 = 102- -N F , (3 2 = -l 2857-—N F + — Nf 

101 5 1249 277+ 180C(3) AT 35 2 
71 = 12-18^' 72 = ^ 108 " 648^' (L?) 

and C(3) = 1.2020569 ... is the Riemann zeta-function. Our initial condition for 
the strong coupling constant will be a s (mz)- Then we will run a s from mz to 
lower scales, i.e. for instance (1q — mz or the upper threshold. On the other 
side, we will run the masses from low to higher scales. We need the inverted 
version of (|1.3|) 



(1.8) 



The beta and gamma functions depend on the number of flavours N p. T here- 
fore, we have to decide whether we have five or four flavours. The trick [20, ^1] [ 



is to built below the heavy quark threshold an effective theory where the heavy 
quark has been integrated out. Imposing agreement of both theories, the full and 
the effective one, at low energies they wrote //. dependent matching conditions 
that express the parameters of the effective theory, with N — 1 quark flavours, 
as a perturbative expansion in terms of the parameters of the full theory with 
N flavours 



Table 3. Running at the NLO and NNLO of the top quark 
mass to m z , a { s 5) {m z ) = 0.118 ± 0.006, af ] {m z ) = 0.117 ± 
0.006. 





fh t {M t ){GeV) 


m t {m t )(GeV) 


m t {m z ){GeV) 


NLO 


167. ±6. 


168. ±6. 


176. ±6. 


NNLO 


165. ±6. 


166. ±6. 


174. ±7. 



1 See also [121]. The mass independent constant coefficient of the < 
dition for the strong coupling constant has been recently corrected by [E 
numerical effect of this correction is very small. 



matching con- 
J|. Nevertheless, the 
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Table 4. Running at the NLO of the bottom quark mass to 
mz and running masses at the running mass scale needed for 
thresholds. For masses extracted from QCD SR a{ 5 \mz) = 
0.118 ±0.006, for lattice ai 5 \m z ) = 0.115 ±0.002 





m b (m b )(GeV) 


m b (m z )(GeV) 


Narison 


4.29 ±0.04 


2.97 ±0.13 


Narison 


4.35 ±0.05 


3.03 ±0.13 


Titard / Yndurain 


4.397±8:S| 


3.07±0.11 


Neubert 


4.37 ±0.09 


3.04 ±0.17 


Crisafulli et al. 


4.17 ±0.06 


2.93 ±0.08 


Gimenez et al. 


4.15 ±0.20 


2.91 ±0.19 


Davies et al. 


4.13 ±0.11 


2.89 ±0.12 


mean 


4.33 ±0.06 


3.00 ±0.12 



Table 5. Running at the NLO of the charm quark mass to 
mz and running masses at the running mass scale needed for 
thresholds. For masses extracted from QCD SR a{ 5 \mz) = 
0.118 ± 0.006, for lattice ai 5) (m z ) = 0.115 ± 0.002 





m c (m c )(GeV) 


m c (m z )(GeV) 


Narison 


1.28 ±0.04 


0.52 ±0.09 


Narison 


1.31 ±0.06 


0.54 ±0.10 


Titard / Yndurain 


l.OUD_ 035 


0.52 ±0.08 


Neubert 


1.17 ±0.12 


0.45 ±0.14 


El-Khadra et al. 


1.36 ±0.19 


0.61 ±0.15 


mean 


1.30 ±0.08 


0.52 ±0.10 



.N-l 



x a 1 ? (a) 

6 7T 



12 ( - + 



llx 



11 



mf( M ) 



(M) 



1 

12 



ox 

y 



89 
36 



(1.9) 



with x = logm 2 (^)//i 2 , where m(/x) is the heavy quark mass which decouples 
at the energy scale fj, and m/(/i) are the lighter quark masses. These matching 
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Table 6. Running of the light quark masses to mz- For masses 
extracted from QCD SR af\m z ) = 0.118 ± 0.006, for lattice* 
ots(mz) = 0.115 ± 0.002. First box is NLO running, second 
and third boxes are NNLO running. 





m s (lGeV)(MeV) m s (m z )(MeV) 


Narison 
Allton et al* 


222. ± 22. 105. ±28. 
156. ±17. 78. ±15. 





m s (lGeV)(MeV) 


m 8 (m z )(MeV) 


Narison 


197. ±29. 


88. ±31. 


Jamin / Miinz 


189. ±32. 


85. ±32. 


Chetyrkin et al. 


171. ± 15. 


75. ± 23. 


mean 


186. ±30. 


83. ±30. 



Bijnens et al. 


Ioffe et al. 


(fh u + m d )(lGeV) = (12.0 ± 2.5)MeV {m d - 
(fh u + m d ){m, z ) = (5.4 ±2. 2)MeV (fh d - 


m u )(lGeV) = (3.±l.)MeV 
m u )(m z ) = {lA±Q.7)MeV 



conditions make the strong coupling constant and the light quark masses discon- 
tinuous at thresholds. However, taking the matching in this way we ensure, as 



pointed out explicitly in [121], that the final result is independent of the particu- 
lar matching point we choose for passing the threshold. As it is independent, the 
easiest way to implement a heavy quark decoupling is to take the threshold scale 
as the running mass at the running mass scale, i.e., nth = rh(fh) or equivalently 
x = 0, then the discontinuity appears only at two loops matching. 

We have summarized in tables [| |], || and [| the results for the quark masses 
running until m z . By NLO we mean connection between the perturbative pole 
mass and the running mass dropping the 0(a 2 s ) term, running to two loops 
and matching at one loop, i.e., strong gauge coupling and masses cont inu ous at 
/ith — m(m). Three loop running and matching as expressed in eq. (|l.9|) with 
x = correspond to NNLO. For consistency with the original works we perform 
the bottom and the charm quark mass running just at NLO. For the light quarks 
the running is consistent to NNLO using the threshold masses, m(m), of the 
bottom and charm quarks determined at NLO. 

We propagate the errors in the running in such a way we maximize them. 
The relative uncertainty in the strong coupling constant decreases in the running 
from low to high energies as the ratio of the strong coupling constants at both 
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scales. On the contrary, in the way we propagate the errors, induced by the 
strong coupling constant error, the quark mass relative uncertainty increases 
following 



i - i w t - i w 270 / ol s {ii) , 

e r {m[m z )) e r {m(fi)) + — — - 1 e r {a s {m z )). 

Pa \a s (m z ) 



(1.10) 



The absolute quark mass propagated error value at high energies depends on 
the balance between the quark mass at low energies and the strong coupling 
constant errors. Case the quark mass error contribution dominates the relative 
uncertainty remains almost fixed and since the running mass decreases at high 
energies its absolute error decreases too. Case the strong coupling constant error 
contribution is the biggest the absolute quark mass uncertainty increases. In this 
last situation, a possible evaluation of the bottom quark mass at the m z scale 
would be considered even competitive to low energy QCD Sum Rules and lattice 
calculations with smaller errors. 

We have depicted in figure |l|. the running of the bottom quark mass, we 
run fh b (fh b ) = 4.39(GeV) with a s (m z ) = 0.112 and m b (m b ) = 4.27(GeV) with 
a s (m z ) = 0.124, the extreme quark mass and strong coupling constant values 
for fh b {fh b ) = (4.33 ± 0.06)GeV and a s {m z ) = 0.118 ± 0.006, and we take 
the difference as the propagated error. In this case, since the strong coupling 
constant uncertainty is the biggest, the absolute bottom quark mass uncertainty 
increases in the running to m z . 

It is informative to notice that the running mass of the top quark is shifted 
about 7(GeV) down from its perturbative pole mass. This shift is of the order 
of its experimental error. Therefore it is important to clarify which mass is 
measured at the CDF and D0 experiments. We have decoupled the top quark 
at m z otherwise it makes no sense to run the top down. This fact shifts down 
slightly the strong coupling constant in m z , from as(m z ) = 0.118 ± 0.006 we 
get a^\m z ) = 0.117 ± 0.006 but has no effect on the masses because the errors 
screen the difference between the theory with 5 and 6 flavours. Curiously the 
running of the top until m z cancels the difference between the running and the 
pole mass, fht(m z ) ~ M t . 

One has to be very careful in comparing the running of the masses obtained 
from QCD Sum Rules and those obtained from lattice because we took differ- 
ent values for the strong coupling constant at m z . Furthermore, we have to 
remember that the error in the running is dominated by the error in the strong 
coupling constant. However it is impressive to notice the good agreement of the 



results obtained in lattice 42, 43] with the running of the masses from QCD 



Sum Rules. Nevertheless, we have to keep in mind that a recent lattice evalua- 



tion [68, 105 1 has enlarged the initial estimated error on the bottom quark mass 
up to 2QQMeV , fh b (fh b ) = (4.15 ± 0.20)GeV, due to unknown higher orders in 
the perturbative matching of the Heavy Quark Effective Theory (HQET) to the 
full theory. 
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Figure 1. Running of the bottom quark mass from low en- 
ergies to the scale. Upper line is the run of rhb(rhb) = 
A.39(GeV) with a s (mz) = 0.112. Bottom line is the run of 
fnb{fhb) = 4.27(GeV r ) with a s (mz) — 0.124. Second picture is 
the difference of both, our estimate for the propagated error. 



We can now play the game of combining the light quark masses of table g 
with the ratios obtained from xPT. The mean value of the strange quark mass 



together with the Bijncns et al. 24 1 result gives 



2TOs =33. ±12., (1.11) 

m u + m d 

in agreement with the xPT result. Being conservative we can also get for the up 
and down quarks m u (lGeV) = (3. ± 2.)MeV and fh d (lGeV) = (9. ± 2.)MeV 
that translate into fh u (m z ) = (1.5 ± l.2)MeV and m d (m z ) = (4.1 ± 1.7)MeV. 

To summarize, the running of the quark masses to the mz energy scale gives 
a running top mass that is around its perturbative pole mass, fh t (mz) = (176. ± 
6.)GeV, for the bottom and the charm quarks we get fhb(mz) = (3.00±0.12)Gey 
and fh c (mz) = (0.52 ± 0.10)GeV respectively, while for the strange quark we 
have a result affected by a big error m s (mz) — (83. ± 30.)MeV. The same 
happens for the up and down quarks, we get rh u (mz) = (1-5 ± 1.2)MeV and 
m d (m z ) = (4.1 ± 1.7)MeV. 
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Three-jet observables at leading order 

As wc mentioned in the previous chapter in the Standard Model of elec- 
troweak interactions all fermion masses are free parameters, and their origin, 
although linked to the spontaneous symmetry breaking mechanism, remains se- 
cret. Masses of charged leptons are well measured experimentally and neutrino 
masses, if they exist, are also bounded. In the case of quarks the situation is 
more complicated because free quarks are not observed in nature. Therefore, one 
can only get some indirect information on the values of the quark masses. For 
light quarks (m q < 1 GcV, the scale at which QCD interactions become strong), 
that is, for u-,d- and s-quarks, one can define the quark masses as the parame- 
ters of the Lagrangian that break explicitly the chiral symmetry of the massless 
QCD Lagrangian. Then, these masses can be extracted from a careful analysis of 
meson spectra and meson decay constants. For heavy quarks (c- and 6-quarks) 
one can obtain the quark masses from the known spectra of the hadronic bound 
states by using, e.g., QCD sum rules or lattice calculations. However, since the 
strong gauge coupling constant is still large at the scale of heavy quark masses, 
these calculations are plagued by uncertainties and nonperturbative effects. 

It would be very interesting to have some experimental information on the 
quark masses obtained at much larger scales where a perturbative quark mass 
definition can be used and, presumably, non-perturbative effects are negligible. 
The measurements at LEP will combine this requirement with very high exper- 
imental statistics. 

The effects of quark masses can be neglected for many observables in LEP 
studies, as usually quark masses appear in the ratio m 2 q /m 2 z . For the bottom 
quark, the heaviest quark produced at LEP, and taking a 6-quark mass of about 
5 GeV this ratio is 0.003, even if the coefficient in front is 10 we get a correction 
of about 3%. Effects of this order are measurable at LEP, however, as we will 
see later, in many cases the actual mass that should be used in the calculations 
is the running mass of the &-quark computed at the mz scale: rhb(mz) ~ 3 GeV 
rendering the effect below the LEP precision for most of the observables. 

While this argument is correct for total cross sections for production of b- 
quarks it is not completely true for quantities that depend on other variables. 
In particular it is not true for jet cross sections which depend on a new variable, 
y c (the jet-resolution parameter that defines the jet multiplicity) and which in- 
troduces a new scale in the analysis, E c = mzy/y^. Then, for small values of 
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y c there could be contributions coming like m 2 /E 2 = (mb/mz) 2 /y c which could 
enhance the mass effect considerably. In addition mass effects could also be en- 
hanced by logarithms of the mass. For instance, the ratio of the phase space for 
two massive quarks and a gluon to the phase space for three massless particles 
is 1 + 8(m q /mz) 2 log(m q /mz)- This represents a 7% effect for m q = 5 GeV and 
a 3% effect for ni q = 3 GeV. 

The high precision achieved at LEP makes these effects relevant. In fact, 
they have to be taken into account in the test of the flavour independence of 
&s( m z) @5 ^ 0> ^ > 41 1 . In particular it has been shown 130| that the biggest 
systematic error in the measurement of a b s (mz) (a 8 obtained from 66-production 
at LEP from the ratio of three to two jets) comes from the uncertainties in the 
estimate of the quark mass effects. This in turn means that mass effects have 
already been seen. Now one can reverse the question and ask about the possibility 
of measuring the mass of the bottom quark, mi,, at LEP by assuming the flavour 
universality of the strong interactions. 

Such a measurement will also allow to check the running of fhb(lJ-) from 
/i = nib to /i = mz as has been done before for a s (n). In addition fhb(mz) is the 
crucial input parameter in the analysis of the unification of Yukawa couplings 
predicted by many grand unified theories and which has attracted much attention 
in the last years 



94 . 



The importance of quark mass effects in Z-boson decays has already been 
discussed in the literature |89, 116]. The complete order a s results for the 
inclusive decay rate of Z — > bb + bbg + ■ ■ ■ can be found]] in [35, 117, 124 1. The 
total cross section and forward-backward asymmetry for e + e~ — > "f*,Z — > qq 
with one- loop QCD corrections were calculated in plf . Analytical results with 
cuts for these quantities were obtained in jnj. The leading quark mass effects 



for the inclusive Z- width are known to order for the vector part 38 ] and 
to order a 2 for the axial- vector part p3|. Quark mass effects for three-jet final 
states in the process e + e~ annihilation into qqg were considered in 
for the photonic channel and extended later to the Z channel in 



79, 


88 


,93 


[112, 


118 



Polarization effects have been studied for instance in [74, 128, 129] for massive 
quarks. In [71 the possibility of the measurement of the quark mass from 



the angular distribution of heavy-quark jets in e + e -annihilation was discussed. 
Recently [13, 



14 



calculations of the three-jet event rates, including mass effects, 
were done for the most popular jet clustering algorithms using the Monte Carlo 
approach. 

We will discuss the possibility of measuring the 6-quark mass at LEP, in 
particular, we will study bottom quark mass effects in Z decays into two and 
three jets. 



1 The order a a corrections to the vector part , including the complete mass dependences, 
were already known from QED calculations [123 . 
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1. The inclusive decay rate Z 



bb 



To calculate at order a s the total decay rate of the Z-boson into massive 
quarks one has to sum up the virtual one-loop gluonic corrections to Z — > bb 
with the real gluon bremsstrahlung, see figure |l|. In addition to renormalized 
UV divergences ^, 1R singularities, either collinear or soft, appear because of 
the presence of massless particles like gluons. Therefore, some regularization 
method for the IR divergences is needed. Bloch-Nordsiek [26| and Kinoshita- 



Lee-Nauenberg [ 82 , 97 theorems assure IR divergences cancel for inclusive cross 



sections. Technically this means, if we use Dimensional Regularization to regu- 



larize the IR divergences [64, 102, 103 1 of the loop diagrams we should express 
the phase space for the tree-level diagrams in arbitrary dimension D = 4 — 2e. 
At order a s and for massive quarks all IR divergences appear as simple poles 
1/e. The IR singularities cancel when we integrate over the full phase space. 



At the end, we obtain the well-known result 47 



m z - 



64tt 



9v 



7T 



-A 



+ 12r 6 ) 

Q/ 

l-6r b + — (1-6^(2 log rj + 1)) 

7T 



(2.1) 



where = mj/m| and gy (gA) are the vector (axial- vector) neutral current 
couplings of the quarks in the Standard Model. At tree level and for the 6-quark 
we have 



9v = -1 



9 A = 1 



(2.2) 



We denote by cw and sw the cosine and the sine of the weak mixing angle. Here 
and below we will conventionally use a s = a s (mz) to designate the value of the 
running strong coupling constant at the ra^-scale. 

It is interesting to note the presence of the large logarithm, log (m^/m|), 
proportional to the quark mass in the axial part of the QCD corrected width, 
eq. (2.1). The mass that appears in all above calculations should be interpreted 
as the perturbative pole mass of the quark. But in principle the expression (2.1) 
could also be written in terms of the so-called running quark mass at the mz 
scale by using 



ml 



m\(mz) 



1 



a 5 A 
2— - 

7T V 3 



ml 



m 



(2.3) 



2 Note that conserved currents or partially conserved currents as the vector and axial 
currents do not get renormalized. Therefore, all UV divergences cancel when one sums properly 
self-energy and vertex diagrams. The remaining poles in e correspond to IR divergences. One 
can see this by separating carefully the poles corresponding to UV divergences from the poles 
corresponding to IR divergences. 
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FIGURE 1. Feynman diagrams contributing to the decay rates 
Z — ► bb, Z — > 66g at order a s . 



Then, we see that all large logarithms are absorbed in the running of the 
mass from the mj, scale to the mz scale pjjfl and we have 

^2 



r 9 
1 6 = m z 2 — 



^(l + ^(l + 12f fc )) 



+si i - 6r-, 



— (i-22f 6 ; 

7T 



where r;, = m^rnxj/mj 
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This result means that the bulk of the QCD corrections depending on the 
mass could be accounted for by using tree-level expressions for the decay width 
but interpreting the quark mass as the running mass. The same point has been 
stressed in [ 106| for the hadronic width of the charged Higgs boson. On the other 
hand, since fhb(mz) ~ 3 GeV is much smaller than the pole mass, « 5 GeV, 
it is clear that the quark mass corrections are much smaller than expected from 
the naive use of the tree-level result with to,;, ss 5 GeV, which would give mass 
corrections at the 1.8% level while in fact, once QCD corrections are taken into 
account, the mass corrections are only at the 0.7% level. 

In fact, this remarkable nice feature of the running mass holds to all orders 
in perturbation theory [37, 39fl . i.e., all the potentially dangerous terms of the 
type to 2 log m 2 /s can be absorbed in the MS scheme. 



2. Jet clustering algorithms 

According to our current understanding of the strong interactions, coloured 
partons, produced in hard processes, are hadronized and, at experiment, one 
only observes colourless particles. It is known empirically that, in high energy 
collision, final particles group in several clusters by forming energetic jets, which 
are related to the primordial partons. Thus, in order to compare theoretical 
predictions with experiments, it is necessary to define precisely what is a jet in 
both, parton level calculations and experimental measurements. 

At order a s , the decay widths of Z into both two and three partons are IR 
divergent. The two-parton decay rate is divergent due to the massless gluons 
running in the loops. The Z-boson decay width into three-partons has an IR 
divergence because massless gluons could be radiated with zero energy. The sum, 
however, is IR finite. Then it is clear that at the parton-level one can define an 
IR finite two-jet decay rate, by summing the two-parton decay rate and the IR 
divergent part of the three-parton decay width, e.g. integrated over the part 
of the phase space which contains soft gluon emission [ 125[ |. The integral over 



the rest of the phase space will give the three-jet decay rate. Thus we need to 
introduce a "resolution parameter" in the theoretical calculations in order to 
define IR-safe observables. Obviously, the resolution parameter, which defines 
the two- and the three-jet parts of the three-parton phase space should be related 
to the one used in the process of building jets from real particles. 

In the last years the most popular definitions of jets are based on the so-called 
jet clustering algorithms. These algorithms can be applied at the parton level 
in the theoretical calculations and also to the bunch of real particles observed 
at experiment. It has been shown that, for some of the algorithms, the passage 
from partons to hadrons (hadronization) does not change much the behaviour 
of the observables ]2g] , thus allowing to compare theoretical predictions with 
experimental results. In what follows we will use the word particles for both 
partons and real particles. 
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Table 1. The jet-clustering algorithms 



Algorithm 


Resolution 


Combination 


EM 


2( Pl - Pj )/s 


Pk — Pi + P] 


JADE 


2{E l E j ){\ - cos0ij)/s 


Pk = Pi + Pj 


E 


(pi+Pj) 2 /s 


Pk=Pi+ Pj 


DURHAM 


2min(Ef,E?)(l - cos%)/s 


Pk=Pi+ Pj 



In the jet-clustering algorithms jets are defined as follows: starting from a 
bunch of particles with momenta pi one computes, for example, a quantity like 



Uij = 2 ElEj (1 - cos% y 



for all pairs (i, j) of particles. Then one takes the minimum of all y^ and if 
it satisfies that it is smaller than a given quantity y c (the resolution parameter, 
y-cut) the two particles which define this are regarded as belonging to the 
same jet, therefore, they are recombined into a new pseudoparticle by defining 
the four-momentum of the pseudoparticle according to some rule, for example 

Pk =Pi+ Pj ■ 

After this first step one has a bunch of pseudoparticles and the algorithm can 
be applied again and again until all the pseudoparticles satisfy j/y > y c . The 
number of pseudoparticles found in the end is the number of jets in the event. 

Of course, with such a jet definition the number of jets found in an event 
and its whole topology will depend on the value of y c . For a given event, larger 
values of y c will result in a smaller number of jets. In theoretical calculations 
one can define cross sections or decay widths into jets as a function of y c , which 
are computed at the parton level, by following exactly the same algorithm. This 
procedure leads automatically to IR finite quantities because one excludes the 
regions of phase space that cause trouble. The success of the jet-clustering 
algorithms is due, mainly, to the fact that the cross sections obtained after 
the hadronization process agree quite well with the cross-sections calculated at 
the parton level when the same clustering algorithm is used in both theoretical 
predictions and experimental analyses. 
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There are different successful jet-clustering algorithms and we refer to 23 , 



for a detailed discussion and comparison of these algorithms in the case of 
massless quarks. 

In the following we will use the four jet-clustering algorithms listed in the 
table [j], where y/s is the total centre of mass energy. In addition to the well- 
known JADE, E and DURHAM algorithms we will use a slight modification 
of the JADE scheme particularly useful for analytical calculations with massive 
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The phase-space for qqg 




Figure 2. The phase space for Z — > b{pi)b{p2)g{pz) in the 
plane y\ = 2(pi -p^)/ s and y-i — 2(pi -pz)l 's with cuts (t/ c = 0.04 
and y c — 0.14) for the different algorithms. The mass of the 
quark has been set to 10 GeV to enhance mass effects in the 
plot. 



quarks [25 



It is defined by the two following equations 
_ 2{ Pi ■ Pj ) 

yn - s 

and 

Pk =Pi+Pj- 

We will denote this algorithm as the EM scheme. For massless particles and at 
the lowest order E, JADE and EM give the same answers. However already at 
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order a 2 s they give different answers since after the first recombination the pseu- 
doparticles are not massless anymore and the resolution functions are different. 

For massive quarks the three algorithms, E, JADE and EM are already dif- 
ferent at order a s . The DURHAM (Kt) algorithm, which has been recently 
considered in order to avoid exponentiation problems present in the JADE al- 
gorithm [^3|, EjlL ^] , is of course completely different from the other algorithms 
we use, both in the massive and the massless cases. 

In figure || we plotted the phase-space for two values of y c (y c — 0.04 and 
y c = 0.14) for all four schemes (the solid line defines the whole phase space for 
Z — > qqg with m q = 10 GeV). 

There is an ongoing discussion on which is the best algorithm for jet clus- 
tering in the case of massless quarks. The main criteria followed to choose them 
are based in two requirements: 

1. Minimize higher order corrections. 

2. Keep the equivalence between parton and hadronized cross sections. 

To our knowledge no complete comparative study of the jet-clustering algorithms 
has been done for the case of massive quarks. The properties of the different 
algorithms with respect to the above criteria can be quite different in the case 
of massive quarks from those in the massless case. The first one because the 
leading terms containing double- logarithms of y-cut (log 2 (y c )) that appear in 
the massless calculation (at order a s ) and somehow determine the size of higher 
order corrections are softened in the case of massive quarks by single-logarithms 
of y c times a logarithm of the quark mass. The second one because hadroniza- 
tion corrections for massive quarks could be different from the ones for massless 
quarks. 

Therefore, we will not stick to any particular algorithm but rather present 
results and compare them for all the four algorithms listed in the table [l]. 

3. Two- and three-jet event rates 

At the parton level the two-jet region in the decay Z — > b{jp\)b{p2)g{p3) is 
given, in terms of the variables [j 2/13, 2/23 and 2/12, by the following conditions: 

2/13 < Vc or 2/23 < Vc or y 12 < y c . (2.5) 

This region contains the IR singularity, 2/13 = 2/23 — and the rate obtained 
by the integration of the amplitude over this part of the phase space should be 
added to the one-loop corrected decay width for Z — ► bb. The sum of these two 
quantities is of course IR finite and it is the so-called two-jet decay width at 
order a s . The integration over the rest of the phase space defines the three-jet 
decay width at the leading order. It is obvious that the sum of the two-jet and 
three-jet decay widths is independent of the resolution parameter y c , IR finite 



3 Only two are independent, in the EM algorithm j/12 = 1 — 2r\, — yri — J/23- 
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and given by the quantity T b = T(Z — > bb + bbg + ■ ■ ■ ) calculated in section [l]. 
Therefore we have 

r b =T b 2j (y c )+T b 3j (y c ) + -- - . (2.6) 

Clearly, at order a s , knowing T b and T%j(y c ) we can obtain r|,(j/ c ) as well. 

The calculation of Tg -(y c ) at order a s is a tree-level calculation and does not 
have any IR problem since the soft gluon region has been excluded from phase 
space. Therefore the calculation can be done in four dimensions without trouble. 
The final result can be written in the following form 

r 3, = ™^^p^ V {9v H v\vc,r b ) +g 2 A Hf{y c ,r b )) , (2.7) 

where the superscript (°) in the functions Hy, A s{y e ,Tb) reminds us that this is 



only the lowest order result. Obviously, the general form 2.7 is independent of 
what particular jet-clustering algorithm has been used. 

In the limit of zero masses, r b — 0, chirality is conserved and the two func- 
tions Hy\y c ,r b ) and H A °\y c ,rb) become identical 

H$\y c ,0) = H^(y c ,0) ee A^(y c ) . (2.8) 

In this case we obtain the known result for the JADE-type algorithms, which is 
expressed in terms of the function (y c ) given by []. 



6y c - -y c + (3 - 6y c ) log 



3 2" 2" c ' °\l-2y 



21og 2 (-^)+4lV Vc 



, / vi (2 ' 9) 

The A(°\y c ) function is also known analytically for the DURHAM algorithm 
lj j33| . Analytical expressions for the functions Hy\y c ,r b ) and H A °\y c ,r b ) 



are given for the EM algorithm in 25 



To see more clearly the size of mass effects we arc going to study the following 
ratio of jet fractions 

( H ( v °\y c ,r b ) H A °\y c ,r b ) \ ( a . 



where we have defined 



(2.10) 



CV - ^ C A - — ^ 

V ,2 _i_ „2 ' ^ A „2 i „2 



ffv +9a 9v+ 9a 

In eq. ( [2.10| ) we have kept only the lowest order terms in a s and r b . The last 
factor is due to the normalization to total rates. This normalization is important 
from the experimental point of view but also from the theoretical point of view 



'Note that with our normalization A^°'(y c ) = \ A{y c ), with A(y c ) defined in 
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because in these quantities large weak corrections dependent on the top quark 



mass 
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cancel. Note that, for massless quarks, the ratio T d j{y c )/Y d 
is independent on the neutral current couplings of the quarks and, therefore, it 
is the same for up- and down-quarks and given by the function A (0) . This means 
that we could equally use the normalization to any other light quark or to the 
sum of all of them (including also the c-quark if its mass can be neglected) . 



The ratios R 



bd 















DURHAM 







0.05 0.1 0.15 0.2 
Vc 



Figure 3. The ratios R b 3 d (see eq. ( 2.10 )) for the four algo- 
rithms. Solid lines correspond to rn;, = 5 GeV, dashed lines 
correspond to mt = 3 GeV and dotted lines give our estimate 
of higher order corrections to the m& = 5 GeV curve. For the 
JADE algorithm we have also included the results of the analysis 
of the data collected during 1990-1991 by the DELPHI group (§. 
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3.1. Estimate of higher order contributions. All previous results come 
from a tree-level calculation, however, as commented in the introduction, we do 
not know what is the value of the mass we should use in the final results since 
the difference among the pole mass, the running mass at (i = m b or the running 
mass at fi = mz are next-order effects in a s . 

In the case of the inclusive decay rate we have shown that one could account 
(with very good precision) for higher order corrections by using the running 
mass at the mz scale in the lowest order calculations. Numerically the effect of 
running the quark mass from mj, to mz is very important. 

One could also follow a similar approach in the case of jet rates and try 
to account for the next-order corrections by using the running quark mass at 
different scales. We will see below that the dependence of R^ d on the quark 
mass is quite strong (for all clustering schemes); using the different masses (e.g. 
nib or rhb(Mz)) could amount to almost a factor 2 in the mass effect. This 
suggests that higher order corrections could be important. Here, however, the 
situation is quite different, since in the decay rates to jets we have an additional 
scale given by y c , E c = mz^fy~ c , e.g. for y c = 0.01 we have E c — 9 GeV and for 
y c = 0.05, E c = 20 GeV. Perhaps one can absorb large logarithms, log(mb/mz) 
by using the running coupling and the running mass at the /i = mz scale, but 
there will remain logarithms of the resolution parameter, log(y c ). For not very 
small y c one can expect that the tree-level results obtained by using the running 
mass at the mz scale are a good approximation, however, as we already said, 
the situation cannot be settled completely until a next-to-leading calculation 
including mass effects is available. 

Another way to estimate higher order effects in R ^ d is to use the known 
results for the massless case [23, 55,|57, 84 , |86| , 87,|9~1, 122 1. Including higher 



order corrections the general form of eq. ( |2.7| ) is still valid with the change 
H$ A Jy c ,rb) — > H V ( A) {y c ,r b ). Now we can expand the functions H V ( A) (y c ,r b ) 
in a s and factorize the leading dependence on the quark mass as follows 



H v(A) (y c ,r b )=AM(y c ) + ^ A^(y c ) 

+ r b (B { °\ A) (y c ,rb) + ^B$ A) (y c ,rbj) +■■■ . (2.11) 

In this equation we already took into account that for massless quarks vector 
and axial contributions arc idcnticalfl 



5 This is not completely true at 0(a^) because the triangle anomaly: there are one-loop 
triangle diagrams contributing to Z — ► bbg with the top and the bottom quarks running in the 
loop. Since rat ^ rn^ the anomaly cancellation is not complete. These diagrams co ntrib ute to 
the axial part even for m;, = and lead to a deviation from Ay\y c ) = (y c ) |73[ . This 
deviation is, however, small |73] and we are not going to consider its effect here. 
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Then, we can rewrite the ratio R^ d , at order a s 



R b 3 d = 



1 + n I Cy 



A«% c ) 



1 + ^ 

7T 



as follows 



B$\y c ,r b ) AM(y. 



B$\y c ,r b ) AW(y, 



-CA- 



Bf{y c ,n) 
A(°>(|fc) 



sW(y c ,^) A«(y c ) 



6r b |ca(1 



a (2/c,r 6 ) 




7T 



A(°)(y c ) 
logr 6 ) - cy2 — } 

7T J 



(2-12) 



The functions By\y c ,r b ) and B A '(y Cl rb) where calculated in [|25| . The 



,(0) 



lowest order function for the massless case, A^°\y c ), is also known analytically 
for JADE- type algorithms, eq. fl2.9|) and refs. [ 23| , 



911, and for the DURHAM 



algorithm pi, 33 . A parametrization of the function A^(y c ) can be found 
in [23] for the different algorithms [| As we already mentioned this function 
is different for different clustering algorithms. The only unknown functions in 
eq. ( |2. 12 ) are By^ (y c , rb) and By (y c ,rb), which must be obtained from a com- 
plete calculation at order a 2 s including mass effects (at least at leading order in 
rb). The second part of this thesis is devoted to this calculation. 

Nevertheless, in order to estimate the impact of higher order corrections in 
our calculation we will assume for the moment that By \{y c , r b)/By A (y c , rb) -C 
AW(y c )/AW(y c ) and take A^(y c )/A^(y c ) fronf] §|,||]. Of course this does 
not need to be the case but at least it gives an idea of the size of higher order 
corrections. We will illustrate the numerical effect of these corrections for R\ d in 
the next subsection. As we will see, the estimated effect of next-order corrections 
is quite large, therefore in order to obtain the &-quark mass from these ratios the 
calculation of the functions By A (y c ,rb) is mandatory. 



3.2. Numerical results for i?§ for different clustering algorithms. 

To complete this section we present the numerical results for R^ d calculated 
with the different jet-clustering algorithms. For the JADE, E and DURHAM 
algorithms we obtained the three-jet rate by a numerical integration over the 
phase-space given by the cuts (see fig. ||). For the EM scheme we used our 
analytical results 25 1 which were also employed to cross check the numerical 
procedure. 

In fig. U we present the ratio R^ d , obtained by using the tree-level expression, 
eq. ( [2.10 ), against y c for nib — 5 GeV and nib — 3 GeV. We also plot the results 
given by eq. Q2.12| ) (with B^ A (y c ,r b )/B^ A (y c ,r b ) = 0) for m b = 5 GeV, which 
gives an estimate of higher order corrections. For y c < 0.01 we do not expect 
the perturbative calculation to be valid. 



6 With our choice of the normalization A^(y c ) = B(y c )/4, where B(y c ) is defined in [p3|] . 
7 For the EM algorithm this function has not yet been computed. To make an estimate of 
higher order corrections we will use in this case the results for the E algorithm. 
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Table 2. Results of the three parameter fits of the functions 
B^(y c ,r b )/A(°)(y c ) = £Lo4"i log"y c in the range 0.01 < 
Vc < 0.20 



Algorithm 


1.(0) 

Ky 




k {2) 

Ky 


1.(0) 
K A 


A 


A 


EM 


-28.58 


-14.64 


-2.72 


-30.67 


— 13.54 


— 2.61 


JADE 


-13.25 


-5.19 


-2.01 


-15.42 


— 4.13 


— 1.90 


E 


25.97 


19.04 


4.68 


23.39 


19.81 


4.71 


DURHAM 


-12.70 


-4.76 


-1.69 


-15.48 


-4.28 


-1.65 



Table 3. Five parameter fits to the LO massless function 
A {0) {y c ) = £Lo fc(n) log"y c in the range 0.01 < y c < 0.20 



Algorithm 








k& 




EM/JADE/E 


1.19248 


2.64729 


1.359065 


-.0173022 


-.001971025 


DURHAM 


.8625 


1.835 


1.0745 


.133 


.01255 



As we see from the figure, the behaviour of R^ d is quite different in the 
different schemes. The mass effect has a negative sign for all schemes except 
for the E-algorithm. For y c > 0.05 the mass effects are at the 4% level for 
nit, — 5 GeV and at the 2% level for nib = 3 GeV (when the tree level expression 
is used). Our estimate of higher order effects, with the inclusion of the next- 
order effects in a s for massless quarks, shifts the curve for mb = 5 GeV in 
the direction of the 3 GeV result and amounts to about of 20% to 40% of the 
difference between the tree-level calculations with the two different masses. For 
both E and EM schemes we used the higher order results for the E scheme. 

To simplify the use of our results we present simple fits in \ogy c to the 
ratios By\(y c ,rb)/A^(y c ), which define R^ d at lowest order, for the different 
clustering algorithms. We use the following parametrization: 

B^ A (y c ,r b )/A^(y c ) = £ fc£ A log"y c . (2.13) 

71=0 

The results of the fits for the range 0.01 < y c < 0.20 are presented in table ||. 

For completeness we include also in table /reftableAO a five parameter fit 
to the LO massless function A^(y c ). For the DURHAM algorithm we take the 



result from [23|. For the EM, JADE and E algorithms that give the same answer 



we have fitted the analytical expression of Eq. (2.E) 
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The functions B^/A^ and B^/A^ 




Figure 4. The functions B { y ] /A^ and B { ° } /A^ for the four 
algorithms. Dashed lines for mt = 3 GeV, dotted lines for 
m,b = 5 GeV and solid lines for our three-parameter fit, 
eq. 



In fig. || we plot the ratios By A {y c , r^,) / A^ (y c ) as a function of y c for the 
different algorithms (dashed lines for = 5 GeV, dotted lines for mi, — 3 GeV 
and solid curves for the result of our fits). As we see from the figure the remnant 
mass dependence in these ratios (in the range of masses we are interested in 
and in the range of y c we have considered) is rather small and for actual fits we 
used the average of the ratios for the two different masses. We see from these 
figures that such a simple three-parameter fit works reasonably well for all the 
algorithms. 

Concluding this section we would like to make the following remark. In 
this chapter we discuss the -Z-boson decay. In LEP experiments one studies 
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the process e + e~ — > (Z7*) — > 66 and, apart from the resonant Z-exchange 
cross section, there are contributions from the pure 7-exchange and from the 
7 — Z-interference. The non-resonant 7-exchange contribution at the peak is 
less than 1% for muon production and in the case of 6-quark production there 
is an additional suppression factor Q\ = 1/9. In the vicinity of the Z-peak the 
interference is also suppressed because it is proportional to Qb{s~m 2 z ) (y/s is the 
e + e~ centre of mass energy). We will neglect these terms as they give negligible 
contributions compared with the uncertainties in higher order QCD corrections 
to the quantities we are considering. 

Obviously, QED initial-state radiation should be taken into account in the 
real analysis; the cross section for 6-pair production at the Z resonance can be 
written as 

T65(«) = / ^i(s')F(s'/ S )ds' (2.14) 

where F(s'/s) is the well-known QED radiator for the total cross section [ jlT| ] 
and, the Born cross section, neglecting pure 7 exchange contribution and the 
7 — Z-interference, has the form 

UttTJ\ s 

a ^ s >- m | ( s _ m |)2 +TO | r | v- 10 ) 

with obvious notation. Note that in this expression can be an inclusive 
width as well as some more exclusive quantity, which takes into account some 
kinematical restrictions on the final state. 



4. Discussion and conclusions 

In this chapter we have presented a theoretical study of quark-mass effects 
in the decay of the Z-boson into bottom quarks at LO in the strong coupling 
constant. Furthermore, we have analyzed some three-jet observables which are 
very sensitive to the value of the quark masses. 

For a slight modification of the JADE algorithm (the EM algorithm) we 



were able to calculate analytically [25| the three-jet decay width of the Z-boson 
into 6-quarks as a function of the jet resolution parameter, y c , and the 6-quark 
mass. The answer is rather involved, but can be expressed in terms of elementary 
functions. Apart from the fact that these analytical calculations are interesting 
by themselves, they can also be used to test Monte Carlo simulations. For the 
EM, JADE, E and DURHAM clustering algorithms we have obtained the three- 
jet decay width by a simple two-dimensional numerical integration. Numerical 
and analytical results have been compared in the case of the EM scheme. 
We discussed quark-mass effects by considering the quantity 
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which has many advantages from both the theoretical and the experimental 
point of views. In particular, at lowest order, the function F(mb,y c ) is almost 
independent on the quark mass (for the small values of the mass in which we are 
interested in) and has absolute values ranging from 10 to 35 (depending on y c 
and on the algorithm), where the larger values are obtained for y c of about 0.01. 

At the lowest order in a s we do not know what is the exact value of the quark 
mass that should be used in the above equation since the difference between the 
different definitions of the 6-quark mass, the pole mass, mb » 5 GeV, or the 
running mass at the m^-scale, fhb{mz) ~ 3 GeV, is order a s . Therefore, we 
have presented all results for these two values of the mass and have interpreted 
the difference as an estimate of higher order corrections. Conversely one can keep 
the mass fixed and include in F(rrib, y c ) higher order corrections already known 
for the massless case. According to these estimates the 0(a s ) corrections can 
be about 40% of the tree- level mass effect (depending on the clustering scheme) , 
although we cannot exclude even larger corrections. 

By using the lowest order result we find that for moderate values of the 
resolution parameter, y c w 0.05, the mass effect in the ratio R^ d is about 4% 
if the pole mass value of the 6-quark, nib « 5 GeV, is used, and the effect 
decreases to 2% if mb — 3 GeV. However, in order to extract a meaningful value 
of the 6-quark mass from the data it will be necessary to include next-to-leading 
order (NLO) corrections since the leading mass effect we have calculated does 
not distinguish among the different definitions of the quark mass (pole mass, 
running mass at the m& scale or running mass at the mz scale). Next chapters 
are devoted to the calculation of the NLO corrections. 

Finally, if gluon jets can be identified with enough efficiency |104| another 
interesting three-jet observable very sensitive to the bottom quark mass is the 
angular distribution 

§ ~ T b dd I T d d$ ' 
where i? is the minimum of the angles formed between the gluon jet and the quark 



and antiquark jets. The angular distribution F&J- was studied at LO in |25|]. We 



S3 

will leave its analysis at NLO order for a future work 



CHAPTER 3 



Transition amplitudes at next-to-leading order 

In the previous chapter we have seen that some three-jet observablcs, in 
particular the following ratio of three-jet decay rate fractions 



are very sensitive to the quark masses and their study at LEP can provide a 
very interesting experimental information about the bottom quark mass. Among 
other applications, such study would allow to perform the first measurement of 
the bottom quark mass far from threshold and, what is more important, it would 
provide, for the first time, a check of the running of quark masses from scales of 
the order of /i ~ 5(GeV) to fj, = mz in the same way the running of the strong 
coupling constant has been checked before. 

However, as we mentioned, in order to extract a meaningful value of the 
bottom quark mass from LEP data it is necessary to include next-to-leading 
order (NLO) strong corrections in since the LO QCD prediction does not 
allow to distinguish among the possible theoretical definitions of the quark mass 
(pole mass, running mass at the nib scale or running mass at the mz scale) that 
numerically are quite different. 

At the NLO we have contributions from three- and four-parton final states. 
The three-jet cross section is obtained by integrating both contributions in the 
three-jet phase space region defined by the jet clustering algorithms considered 
in the previous chapter: EM, JADE, E and DURHAM. This quantity is infrared 
finite and well defined, however the three- and four-parton transition amplitudes 
independently contain infrared singularities. Therefore, some regularization pro- 
cedure is needed. We use Dimensional Regularization because it preserves the 
QCD Ward identities. Since at this order we have diagrams with a three gluon 
vertex it is not possible anymore to regulate the infrared divergences with a 
gluon mass as it would be possible at the lowest order. 

The three-parton transition amplitudes can be expressed in terms of a few 
scalar one-loop integrals. After UV renormalization we obtain analytical expres- 
sions for the terms proportional to the infrared poles and for the finite contribu- 
tions. The finite contributions are integrated numerically in the three-jet region. 
The infrared poles will cancel against the four-parton contributions. 




i%(y c )/r" 
rl-(2/ c )/r d 



•d • 



(3.1) 
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The four-parton transition amplitudes are splited into a soft part in the 
three-jet region and a hard contribution. The soft terms are integrated ana- 
lytically in arbitrary D dimensions in the region of the phase space containing 
the infrared singularities. We obtain analytical expressions for the infrared be- 
haviour of the four-parton transition amplitudes and we show how these infrared 
contributions cancel exactly the infrared singularities of the three-parton tran- 
sition amplitudes. The hard terms are calculated in D = 4 dimensions. The 
remaining phase space integrations giving rise to finite contributions are per- 
formed numerically. 

In this chapter we consider the NLO corrections to the three-jet decay rate 
of the Z-boson into massive bottom quarks. First we present and classify the 
matrix elements of the virtual corrections to the process 

Z(q)^b( Pl ) + b( P2 )+g( P3 ) , (3.2) 

and the tree-level processes 

Z(q) b(p x ) + 6(pa) + g(p 3 ) + g{Pi) , 

Z(q) -> b( Pl ) + b(p 2 ) + q(p 3 ) + q(Pi) , ^ g . 

Z{q) b(pi) + b{p 2 ) + b(ps) + b{ Pi ) , 



where q stands for a light quark and the symbols in brackets denote the particle 
momenta. In the following we will denote by pij the sum of the momenta of 
particles labelled % and j, p^ = pi +Pj. In next chapter we calculate the singular 
IR pieces of these matrix elements in the three-jet region and show how after 
UV renormalization IR singularities cancel among the one-loop corrected width 



of Z — > bbg and the tree-level processes (3.3) 



We use Dimensional Regularization to regularize both the UV and the IR 



divergences [|64, 102, 103 1 with naive anticommutating 75 prescription, see sec- 
tion U of Appendix [Al. A ll the Dirac algebra is performed in D-dimensions with 
the help of the HIP~^Fsj[ Mathematica 2.0. package. We work in the Feynman 
gauge. 

1. Virtual corrections 

The radiative corrections to the process 

Z(q)^b( Pl ) + b(p 2 )+g(p 3 ) , (3.4) 

are shown in figure |]. They contribute to the three-jet decay rate at 0{a 2 s ) 
through their interference with the lowest order Feynman diagrams TO (a) and 
T0(b). We have not depicted the diagrams with selfenergy insertions in the 
quark and gluon external legs since their contribution is just the wave function 
renormalization constant times the square of the lowest order matrix elements. 

Only the one loop Feynman diagrams Tl, T2, T3 and T4 hold UV diver- 
gences, the one-loop box integrals of diagrams T5 and T6 are UV finite since they 
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FIGURE 1. Radiative corrections to the process Z — > bbg. Self- 
energies in external legs have not been shown. Interference with 
the lowest order Feynman diagrams, T0(a) and T0(b), gives rise 
to the 0{a 2 s ) correction. 

contain four propagators. Diagrams T3 and T4 are responsible for the renormal- 
ization of the strong coupling constant and diagram T2 renormalize the quark 
mass. The UV divergences of diagram Tl are cancelled by the UV piece of the 
quark wave function renormalization constant. We perform the UV renormal- 
ization in a mixed scheme where the strong coupling constant is renormalized in 
the MS scheme while the rest is renormalized in the on-shell scheme. After UV 
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renormalization Tl, T2 and T3 become completely finite, T5 remains with only 
simple IR poles whereas T4 and T6 contain up to double IR divergences. Colour 
can be used to preliminary classify the one-loop diagrams. The interference of 
diagrams Tl and T2 with the lowest order diagrams TO carries a colour factor 
Cfi, where Cf = 4/3. Diagrams T3 and T5 generate a Cf{Cf — Nc/2) colour 
factor, with Nc — 3 the number of colours, whereas T4 and T6 produce Cf^c- 
We found the most efficient way to calculate the loop diagrams is to directly 
perform all the traces over the Dirac matrices on the interference matrix elements 
to reduce them to a Lorentz scalar, before integrating over the virtual gluon 
loop momentum k. After this, we end up with a series of loop integrals with 
scalar products of the internal momentum k and the external momenta in the 
numerator. All these vector and tensor loop integrals can be reduced to simple 
scalar one loop integrals by following the Passarino-Veltman reduction procedure, 
see section || of Appendix |X| At the end, our problem is simplified to the 
calculation of only five one loop three propagator integrals, two infrared divergent 
scalar box integrals and some well known scalar one- and two-point functions, 
see Appendix |A|. Moreover, the infrared divergent piece of the two four-point 
functions can be written in terms of two of the three-point one-loop integrals. 



2. Emission of two real gluons 

The calculation of the transition probability for the process 

Z{q) -> b(jn) + b{ P2 ) + g(p 3 ) + g{ P4 ), (3.5) 

from the eight diagrams shown in figure ^| contains in principle 36 terms. How- 
ever, many of them are related by interchange of momentum labels and, at the 
end, only 13 transition probabilities need to be calculated. We follow the nota- 



tion of [55]. The interference of graph Bi with Bj is written as 5-y. The relevant 
13 transitions probabilities, which we display in fig§, can be classified into three 
different subsets depending on the colour factor. In tabic [l] we give the mo- 
mentum label interchanges necessary to generate all the transitions probabilities 
from the thirteen which we choose to calculate. 

Therefore, it is sufficient to consider the following combinations of transitions 
amplitudes 

Class A = -Bll + 2521 + B22 + 532, 

Class B = -541 + 2542 + B62 + 552, (3.6) 

Class C = 2 (£71 + B72 + 582) + ^(B77 + B87), 

plus the interchanges (1 *-* 2), (3 <-> 4) and (1 «2) (3 «-> 4). 

Since only gluons attached to external legs can generate IR divergences in 
the three-jet region we can see immediately that B32 and B52 are fully finite, 
B21, 522, 542 and 562 are IR only in gluon labelled as 3 while 511 and 541 
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Figure 2. Feynman diagrams contributing to the process Z — -> bbgg. 

Table 1. The interchange table relating the graphs for Z — > 
bbgg. 



label 
permutation 


Class A 
C| 


Class B 
C F (C F -\N C ) 


Class C 
C F N C 


(1^2) 
(3-4) 
(1 <-> 2) (3 <-» 4) 


Bll B21 B22 B32 
B44 B64 B66 B65 
B44 B54 B55 B65 
Bll B31 B33 B32 


B41 B42 B62 B52 
B41 B61 B62 B63 
B41 B51 B53 B52 
B41 B43 B53 B63 


B71 B72 B82 B77 B87 
B84 B86 B76 B88 B87 
B74 B75 B85 B77 B87 
B81 B83 B73 B88 B87 



are IR in both gluons, 3 and 4. All the matrix elements of Class A and B 
are free of quark-gluon collinear divergences since we are working with massive 
quarks. Therefore, we will find there just IR simple poles, 1/e, because only soft 
singularities remain. On the other hand, we expect IR double poles for diagrams 
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of Class C because the gluon-gluon collincar divergences are still preserved at 
the three gluon vertex. This argument is not completely true for the transition 
probabilities B77 and B87. Both diagrams individually contain IR double poles 
but when we take into account all the momenta interchanges, i.e., in the square 
of the sum of diagrams B7 and B8, double poles cancel. As we will see, these 
diagrams are related, by Cutkowsky rules, to the diagram with a selfenergy 
insertion in an external gluon leg. That is the reason why only simple IR poles 
can appear. 




B41 B42 B62 




FIGURE 3. Transitions probabilities for the process Z — > bbgg 
On-shcll particles are indicated by the dashed line and the num- 
bers refer to the momentum labels. All the other transition 
probabilities can be obtained by permutation of momentum la- 
bels. 

We must sum only over the two physical polarizations of the produced gluons. 
This is most easily accomplished by summing over the polarizations with 



pol 



(3.7) 
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but including B7- and B8-like Feynman diagrams with "external" ghosts to take 
into account of the fact that the gluon current is not conserved. 

3. Emission of four quarks 

Lastly, we calculate the matrix elements for the decay width of the Z-boson 
into four quarks. Two processes have to be considered 

Z(q) -> b(pi) + b(p 2 ) + q{p 3 ) + q(Pi), 

Z(q) 6(p0 + 6(pa) + 6(p 3 ) + 6(p4), 

where g stands for a light quark. 



CI 




C2 




Figure 4. Feynman diagrams contributing to the process Z 
bbqq where q stands for a light quark. 



For the first process we restrict to the case where the pair of bottom- 
antibottom quarks is emitted from the primary vertex, figure ^. Similar Feynman 
diagrams, where the heavy quark pair is radiated off a light qq system, are also 
possible. Despite the fact that b quarks are present in this four fermion final 
state, the natural prescription is to assign these events to the light channels. 
Obviously, they must be subtracted experimentally. This should be possible 
since their signature is characterized by a large invariant mass of the light quark 
pair and a small invariant mass of the bottom system. Since four fermion final 
states, bbqq, originate from the interference between qq and bb induced ampli- 
tudes to assign them to a partial decay rate of a particular flavour is evidently 
not possible in a straightforward manner. 
The massless calculation 
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, 55 
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to which we want to compare our re- 
sults was performed by summing over all the allowed flavours. No ambiguity 
appears in this case. The massless QCD prediction is proportional to the the 
sum over all the squared vector and axial- vector couplings 



{s/vi+ffki), 



(3.9) 



although this factor cancels in the three-jet decay rate ratios. Our choice is the 
appropriate to get 1 for the massive over the massless three-jet fraction ratio 
in the limit of massless bottom quarks. Conclusively, there is no way to solve 
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this ambiguity and hence in the case of four fermion final states the theoretical 
analysis should be tailored to the specific experimental cuts. 

The transition probabilities Cll, C22 and C12 can generate infrared diver- 
gences, soft and quark-quark collinear singularities, in the three-jet region due to 
the light quarks. In principle these divergences can manifest as double poles. [] 
Nevertheless, as in the case of the gluon-gluon collinear divergences of the tran- 
sition probabilities B77, B88 and B87 double poles cancel in the sum because 
these transition amplitudes are related, by Cutkowsky rules, to the diagrams with 
a gluon selfenergy insertion in the Born amplitudes TO (a) and TO (6), where only 
simple infrared poles can appear. All the other transition probabilities we will 
treat in this section are infrared finite. 




FIGURE 5. Feynman diagrams contributing to the decay width 
of the Z-boson into four massive bottom quarks. 



For exact massless quarks. It is also possible to regularize these infrared divergences 
with a small quark mass. In this case, infrared divergences are softened into mass singularities 
and lead to large logarithms in the quark mass, log(m q /fi). Infrared gluon divergences can be 
regulated at lowest order by giving a small mass, A, to the gluons. At next-to-leading order 
we would violate gauge invariance at the three gluon vertex. 
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We consider now the emission of four bottom quarks. As in the case of the 
emission of two real gluons, from the eight diagrams shown in figure |] we should 
compute only twelve transition probabilities because many of the, in principle, 
possible 36 terms are related by interchange of momentum labels. 

Table 2. The interchange table relating the graphs for Z — > 
bibb. 



label 
permutation 


Class D 
CfTr 


Class E 
C F (C F - \N C ) 


Class F 

c F 


(1^3) 
(2-4) 
(1 <-> 3) (2 <-> 4) 


Dll D12 D22 
D55 D56 D66 
D77 D78 D88 
D33 D34 D44 


D18 D25 D15 D28 D17 D26 
D45 D16 D15 D46 D35 D62 
D27 D38 D37 D28 D17 D48 
D36 D47 D37 D46 D35 D48 


D13 D23 D24 
D57 D67 D68 
D57 D58 D68 
D13 D14 D24 



It is sufficient to consider the following combinations of transitions ampli- 
tudes 

Class D = Dll + 2D12 + D22 , 



Class E = -(D15 + D28 + Dll ■ 
Class F = D13 + 2D23 + D2A , 



D26 + 2(D18 + D25)) 



(3.10) 



plus the interchanges (1 <-> 3), (2 <-> 4) and (1 <-> 3) (2 «-> 4). Due to Fermi 
statistics there is a relative minus sign for diagrams D5 to D8 that is reflected 
in the transition probabilities of Class E and ensures that each helicity ampli- 
tude vanish when both fermions (antifermions) have identical quantum numbers. 
Furthermore, we will call 

Class G =| CI + C2 | 2 (3.11) 

that evidently carries the same colour factor as the Class D transition probabil- 
ities, C F T R , where T R = 1/2. 



The transition amplitudes of Class F are called in the literature [39], singlet 
contributions because they contain two different fermion loops and hence can be 
splited into two parts by cutting gluon lines only. The first contribution to the 
vectorial part arises at 0{o? a ) as a consequence of the non-abelian generalization 
of Furry's theorem. Singlet contributions to the axial part appear already at 
0{al). 

Lastly, let's comment on the Class F-like case where a bottom quark is 
running in one of the loops and a light quarks is in the other. This would 
produce a term proportional to the product of the vector-axial couplings of both 
quarks, gAbgAq, and represents the interference of diagrams of figure || with 
those we have not considered. Nevertheless, since the vector-axial coupling of 
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FIGURE 6. Transition probabilities for the process Z — > bbbb. 
On-shell particles are indicated by a short cutting line and the 
numbers refer to the momentum labels. All the other transition 
probabilities can be obtained by exchange of momentum labels. 



up and down quarks are different in sign, qa u = —QAd = 1, their contribution 
will cancel when summing over all the light quarks and only the diagram with 
bottom quarks running on both loops will survive. 

This section concludes our classification of the transition probabilities we 
must compute in order to achieve the NLO QCD prediction for the three-jet 
decay rate of the Z-boson into massive quarks. As we have commented, IR 
divergences are the main problem that appears when we try to compute them. 
In the next chapter we will show how to extract from these transition amplitudes 
the divergent pieces and how to keep away the finite parts. We will integrate 
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analytically the divergent contributions over phase space and we will show how 
infrared divergences are cancelled. 



TRANSITION AMPLITUDES AT NEXT-TO-LEADING ORDER. 



CHAPTER 4 



Infrared cancellations 



In the previous chapter we have presented the NLO virtual and real correc- 
tions to the decay width of the Z-boson into three jets. We were able to reduce 
the contribution of the real corrections to the calculation of a few transition prob- 
abilities and the contribution of the virtual corrections to the calculation of a few 
scalar one- loop n-point functions, see Appendix |b|. After UV renormalization of 
the loop diagrams we still encounter a plethora of IR divergences which should 
cancel with the soft and collinear singularities of the tree level diagrams when 
we integrate over the three-jet region of phase space. The theorems of Bloch- 
Nordsiek [26 and Kinoshita-Lee-Nauenberg [82, 97] guaranty such cancellation. 



We have classified our transition probabilities following their colour factors. 
It is clear that the cancellation of IR divergences can only occur inside groups of 
diagrams with the same colour factor. The problem of the IR cancellation can be 
simplified if we find a criteria to split our transition probabilities into different 
groups such that the cancellation of IR singularities can be shown independently. 
The key is to depict them as bubble diagrams as we did in figure [| and figure ^ 
and perform all the possible cuts that could lead to three-jet final states, see 
figure |l|. 

It is not difficult to convince ourselves that transition probabilities of Class 
A can be related only to the insertion of the quark selfenergy in the external legs 
of the lowest order diagram TO. The transition probabilities of Class B have the 
same IR structure as the one loop diagram T5. Besides, B71, B72 and B82 lie 
in the same group as the one loop diagrams T4 and T6. And finally, 577, B87, 
Cll, C12, Dll and D22 are related with the gluon selfenergy insertion diagram 
at TO. 

In the following sections we are going to show analytically how the cancel- 
lation of the IR divergences occurs following the previous classification. 



1. Soft divergences 

The phase space of n + 1 particles can be written as the product of the n- 
body phase space times the integral over the energy and solid angle of the extra 
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particle |132( | . In arbitrary D — 4 — 2e dimensions we have 

dPS{n + 1) = ^J^-i ^ D ~ 3 dfi ■ (4-1) 

Suppose E 3 is the energy of one soft gluon: E 3 < w where w, with w very small, 
is an upper cut on the soft gluon energy. Let's consider 

i 

E°- 3 dE 3 dQ- ^ , (4.2) 

(2pi • p 3 y 

that appears for instance in the four parton transition amplitudes of Class A, 
where p 3 is the momentum of the gluon and p\ is the momentum of the quark. 
After integration over the trivial angles we find 

4T(§ -1) ]„ ' J_, (E, -p,i) 2 

with p x the modulus of the threemomentum. In the cikonal region we can 
suppose E\ almost independent from E 3 

Ek-^tl-yu), (4.4) 

with j/24 = 2(j>2 -p4)/m 2 z . The integral over E 3 can be easily done and gives a 
simple infrared pole in e. Since the x dependent part is completely finite we can 
expand in e before the integral over this variable is performed. The final result 
we get is as follows 

1 f w 1 



2(2^ J Q ^^V^pW 



1 (4tt) £ w°- 4 1 
16^2r(l-e)L>-4£; i 2 -p2 

Let's consider now 



I M-log|^ + 21og2 
Pi Ei+ Pi 



(4.5) 



E 3 - 3 dE 3 dfl— — , (4.6) 

3 (2 Pl ■ p 3 )(2p 2 ■ p 3 ) ' 1 ' 

that appears in the four parton transition amplitudes of Class B, with pi the mo- 
mentum of the antiquark. In this case we can perform a Feynman parametriza- 
tion over the two momenta scalar products 

If 1 1 

(pi • P3)(P2 ■ P3) Jo [P3 ■ {pi + {P2 - Pi)y)] 2 

We get therefore the same integral structure as in Eq. (4.5), integral over y left, 
but with the fourmomentum 

Pa =Pi + {P2 ~Pi)y ■ (4.8) 
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instead of p\ . For the divergent part we get 

1 



i r F D 



3 dE?<m- 



2{2ir) D ~ 1 Jo d (2pi-P3)(2p2-ps) 



D-4 



1 (4tt) £ to 
~ 16tt 2 r(l - e) D - 4 

where 



2 



(?yi2 + 2r 6 )/3i 



logci2 + 0(e) 



(4.9) 



2/2 fl A 4r b l-/3l2 

r fe = m b m z , f3 12 = . I — — , ci 2 = . 

V yi2 + 2r 6 l + /3i2 ( 41Q ) 

It is not easy to get a full analytical expression for the finite part. Nevertheless, 
the integral over the y parameter can be done numerically, using for instance 
a simple GAUSS integration. Notice this last integral has the same divergent 
structure as the scalar one loop C05 function defined in Eq. (A. 17). 
The matrix element of diagram B2 reads 

B2 = i-^-g 2 s (T a T b )u( Pl ) 

W . 1 7^ , 1 ~ Jk(9v + 5.475) j v(p 2 )e* a v (p 3 )er (P4)4(9) . (4-H) 
L V 13 - m b V 134 - rn b J 

where T a = X a /2 are the colour Gcll-Mann matrices, e*^, e* h a and e z stand 
for the polarization vector of the two gluons and the Z-boson and gy (<7a) are 
the vector (axial- vector) neutral current couplings of the quarks in the Standard 
Model. At tree-level and for the bottom quark we have 

9v = -l + ^s 2 w , g A = l- (4.12) 

We denote by cw and sw the cosine and the sine of the weak mixing angle. 
With the help of the equation of motion we can write 

u(pih v (]/i3 + m b ) = u(p 1 )( 1 v ] / 3 + 2p v 1 ) . (4.13) 

In the limit of gluon labelled as 3 soft, p% — > 0, it reduces to 2p v x u(j>\) and B2 
behaves as the Born amplitude TO (a) 

B2 ~ -g s (T a )M T o( a ) £ *> 3 ) . (4.14) 

2/13 

Hence, in this limit we find for the transition amplitude £?22 the following result 

B22~-g 2 s C F ^ \T0(a) | 2 . (4.15) 

Vl3 

Integrating over phase space with the help of Eq. (4.5) we find 

n (A-kY 2nr 2e 

dPS(4)B22 ~ C F -?--±—!— dPS(3)[l + 0(e)] \ T0(a) | 2 . 

4rrr(l-e) e (41g) 
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Same argument can be applied to the transitions amplitudes _B11 and B21 
for which we get 



Bll 



B21 



- £ (\T0(a) | 2 + 
-T0(a)*TQ{b) 



T0(&) 



(4.17) 



For the moment we don't need to consider B32 because it is infrared finite. Since 
the transition amplitude -Bll can be soft in both gluons, labelled 3 and 4, we 
have splited off the contribution of each one by performing a partial fractioning 



1 



1 



viz 



2/24 



1 

viz 



1 

vL 



(4.18) 



The first term should be integrated in the 1-3 system, second term in the equiv- 
alent 2-4 system. Relabelling in both cases the remaining hard gluon as 3 we 
arrive at the result quoted in Eq. (4.17). 

Finally, taking into account all the momenta permutations, we find for the 
cikonal contribution of the diagrams of Class A to the decay width of the Z-boson 
into massive bottom quarks the following result 

(4tt) £ w~ 2t 1 



r„(ClassA) = |c F g r(1 



2m z 



dPS(3) | T0(a) 



■T0(b) | 2 , 

(4.19) 



where 2! is the statistical factor and 
1 



16c^ 



9 2 9 2 s Cf 
8{D - 2 



9 2 v 



T0(a) +T0(6) | 2 

fin 



4r h ) 



2 2 
^13^23 



+9 2 a(D - 2) 
with 



(1 - 4r b ) 



2 2 
VlzV2Z 



2(D-2) 2 
+ 2(D-2 



2/23 



V2Z 
2/13 



4r h) ( — 

2/23 



2/23 
2/13 



2){D-A) 
K4(D-4- 



(4.20) 
+" 4r 6 ) , 

(4.21) 



hp = 2/132/23(1 - 2/13 - 2/23) - n{yiz + 2/23) , 

is the squared lowest order transition amplitude in D-dimensions. 

The wave function renormalization constant of the quark propagator con- 
tains two pieces 



Z-2 = I—Cf—tq 

4-7T 



A uv - log -± 



(3 - to) 



2 ~i 

Am - log -1 



2(2 - r ) 



(4.22) 



where A = 1/e — 7_e + log47r and To is the gauge parameter, To = 1 since we 
work in the Feynman gauge. The first one is the usual UV piece and cancels the 
UV divergences of diagram IT. The second piece comes from the residue of the 
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MS renormalized propagator in the pole. When we include it in the self-energy 
diagrams of external quarks we get a contribution 

- 2^C F 2A IR -^— [ dPS(3) | T0(a) + TQ(b) | 2 , (4.23) 
4-7T 2m z J 

which exactly cancels the IR diver gences of ( 4.19 ). 

As we made with B2 in ( 4.14 ) the matrix element of diagram BA reads 

B4 ~ g s (T a ) (c F - ^f) M r0 (a)er(p 3 ) , (4.24) 

for gluon labelled as 3 soft. The colour factor comes from the following properties 
of the SU (3) Gell-Mann matrices 

[T a ,T b ] =lfabcT C , 

N c ( 4 - 25 ) 
fabcT b T c = l-^-T a . 

Therefore, the transition amplitude -B42 behaves as 

2 ( n N c\ Mpi-Pa) , ™, s |2 



B42 ^gi[C F -^ _^±L | T Q(a) \> , (4.26) 

V 2 / 2/132/23 

in the eikonal region. Including phase space, from Eq. (4.9) we get for its diver- 
gent piece 

dPS(4)B42 ^(c F -^ as 



2 J An T(l - e) 

2w- 2e 



-dPS(3) 



:h " logci2 + 0(e) 



(yi2 + 2r b )Pi2 



T0(a) | 2 . (4.27) 



As before we find for the other transition amplitudes of Class B the following 
behaviour 

£41 - i[T0(a)*T0(6) +T0(b)*T0(a)} , 

B62 ~ -T0(a)*T0(b) , 
e 

where for _B41 we made the following partial fractioning 

1 1/11 



2/132/232/142/24 2/132/23 + 2/142/24 V 2/132/23 2/142/24 

The final answer for this class of diagrams reads 
r„(ClassB) = A^-^ ^ 



(4.28) 



(4.29) 



2! V 2 J 47rr(l -e) 

2£ 1 f Jn *,^ 2/12 



2mz 



f dPS(3)- yi2 logc 12 | T0(o) + T0(b) | 2 . (4.30) 

J (2/12 + 2r b )Pi2 
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Let's consider now the one-loop box diagram T5 

br r ar r b\ 1 U K 



T5(a) = -giu{ Pl ){T b T a T b ) 



l a [(¥ + ]/i)+ m b ]Y[{¥ + 1/13)+ m^figy + g A lhW -1/2) + m b ]-f' 
k 2 [(k+ Pl ) 2 ~ m 2 b }[(k +P13) 2 - rn 2 b ][{k - p 2 ) 2 - mg] ( 4 3-^ 

/ ' g,vl ^) £ ;> 3 )4(?)- 



6W - (1 - TO) 



fc 2 

Again we apply the equation of motion of the quarks 

(4.32) 

[(¥-]/ 2) + m b }^ v( P2 )^(¥l a -2pl)v(jp 2 ) , 

at both sides of the previous equation. Dropping all the fi factors in the nu- 
merator and expanding the third propagator for small loop momentum we arrive 
to 

T5(a) ~ -^(4pi ■p 2 )C05(y 12 )T0(a) , (4.33) 



where C05 is defined in Eq. (A. 14). It is straitforward to see from the solution 



(A. 17) to the one-loop integral C05 how the interference of the one loop ampli- 



tude of Eq. (4.33) with the lowest order amplitude TO cancels exactly the IR 



divergences of the diagrams grouped in Class B, Eq. (4.30). 

2. Collinear divergences 

To show how the cancellation of the infrared divergences for diagrams of 
Class C occurs is quite more difficult. Since we have gluon-gluon collinear di- 
vergences it is not enough to make a cut in the energy of one of the gluons, as 
we made in the previous section, to extract the divergent piece. In this case the 
appropriate variable over which we have to impose a cut is j/34 = 2(p3 • p4)/m 2 z 
as it includes both kind of divergences, soft and collinear, in the limit of j/34 — > 0. 
Looking at the four partons phase space in the so-called "system 3-4" defined in 
section^ of appendix [B| we can notice several things. First, in the limit 7/34 — > 
the function h p that defines the limits of the phase space reduces to the three 



parton phase space h p function ( P3.4|) . Second, the momentum p 3 4 — p 3 + p t 
behaves as the momentum of a pseudo on-shell gluon because p 2 4 — ► in this 
limit. Therefore, it is possible to factorize from the four-body phase space a 
three-body phase space with an effective gluon of momentum p 34 . Only integra- 
tion over the extra variables should be performed to show the cancellation of the 
infrared divergences. 

2.1. Diagrams containing gluon-gluon collinear divergences. Let's 
consider the following phase space integral 

PS(4) , (4.34) 

yi32/34 
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in the called system 3-4 ( [B.15 ), where as usual p\ denotes the momentum of 
the quark and p 3 and P4 are the two gluon momenta. This integral can be soft 
for gluon labelled as 3 and collinear because of the scalar product of the two 
gluon momenta 2/34- In the eikonal region, 2/34 < w with w very small, we can 
decompose the four body-phase space as the product of a three-body phase space 
in terms of variables 2/134 and 2/234 times the integral over 2/34 and the two angular 
variables v and 0' 

PS(4)—L_ = P5(3)(2/i34, 2/234) 

J/132/34 

T^fTT^ r^ 34 fdv{v[l-v))-^ ^ d& sin" 2 ^'^ . (4.35) 
16tH r(l - e) J J N g > J 2/13 

S = 1/2!, the statistical factor. In the 3-4 system the two-momenta invariant 
2/13 can be written in terms of the integration variables as 

2/13 = ^ (2/134 - V2/134 - 4r b 2/34(1 - 2v)j . (4.36) 

The 2/13 factor is independent of the 8' angle so we can get rid of this first 
integral. Second observation, 2/13 contain a piece independent of the v angle and 
another one that is odd under the interchange v <-*• (1 — v). With the help of 
this symmetry 

1/1 1 



2\A + f(v) ' A-f(v)J = A*-f{y)* ' (4 ' 3?) 
we can avoid the use of the square root of Eq. ( 4.36| ) and write our integral as 

PS(4)— = PS(3)(2/i34, 2/234) J 2 Jf } \ (4-38) 
2/132/34 I67H T(l - e) 

l-e,, f 1 ,,,, \ \ — c 2 yi34 



2/34 £ ^34 / dv(v(l - v)) 



lo ~" V "^ "" yf 3 4- (y 2 134 -4r b y M ){l-2v) 2 
After integration over 2/34 we get the first infrared pole 

dv(v(l - v)) £ 2 ^i[l, -e, 1 



e 22/134 Jo ~ ' 2/i34«(l - u) ( 4 39 ) 

Some simple mathematical manipulations |85|, 113 , 115 ] allow us to express 
the hypergeometric function in terms of a simple dilogarithm 

2 F 1 [l,-e,l-e,~]=(-^—) 2 F 1 [-e,-e,l 



b \ a + b J a + b 



b 



a - 



l + e 2 U 2 ( — )+0(e 3 ) 



(4.40) 
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to finally obtain 
— 1 w 



e 2y 



121 



/ dv(v(l - v))- 1 - 2 ' (v(l -v) + A(l - 2v) 2 ) 
Jo 

A(l-2v) 2 



1 + e z Li 2 



v(l - v) + A(l - 2v) 2 



+ 0(6 3 ) 



(4.41) 



with A = rbw/y 2 34 . The contribution of the dilogarithrn function reduces to 
its value at the border of the integration region, i.e., just additional Lz2(l) = 
7r 2 /6 factors to the finite part. Further contributions would be next order in 
e. Therefore, it is enough to analyze the first two factors. This integrand is 
symmetric under the interchange v <-» (1 — v ). We perform the integral only in 
half of the integration region and apply the following change of variables 



u = 4v(l — v) 



(4.42) 



to get 



-1 w~ e A e 
e 2/134 



>4e ^ 

Jo 



duu 



-l-2e 



(l-«) 



-1/2 



which gives again an hypergeometric function 



-1 1 



e 2/134 Vyi34/ 



r(-2 £ ) 

r(l-2e) 



iFi\ 



1 -~A ]U 



1 1 

-2e,--2e,l- — 
'2 4A J 



(4.43) 



(4.44) 



The hypergeometric function is already 1 + 0(e 2 ) but now it is not possible to 
write it in terms of a simple dilogarithrn function. If we are interested just in 
the pole structure we can stop here. Observe that neither the double pole or 
the simple one depend on the eikonal cut w. To get the complete finite part 
further mathematical manipulations must be performed on the hypergeometric 
function. 
Be 

2 Fi[-e,-2e,l-2e,l-^] = 



(4^T ( i-2e ) | f ^^ 2 F 1 [-2e,i-e,l-MA] 



+ 



r(-e) 



r(-2e)r(i-e) 



-e,^,l + e,4A] 



. (4.45) 
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For w enough small, w <C 2/i34/ r f>, we can use the Gauss series |3| to expand the 
hypergeometric functions around A — > 0. The final result we obtain is as follows 

dPS{4)—!— = dPS(3)(y lu , 2/234) 

2/132/34 

S W l/l 1 , n 1 2 r b 1 2 rr 2 
777 + 77 lo § 77" +7 lo S 77- - 7 lo § A - 7T 



16tt 2 T(l - e) y 134 Ue 2 2e * y 2 34 4 b yf 34 2*4 

-2AlogA-yl 2 (2 + 31og^) + 0(A 3 ) 



Notice we get the same infrared poles as in the one-loop three point function 
C03, Eq. (A. 16), appearing in the one-loop amplitudes T4 and T6. 

2.2. Gluon selfenergy-like diagrams. We discuss now the infrared be- 
haviour of the transition amplitudes B77, B87, Cll, C12 and those grouped in 
Class D. Diagrams of Class D, that correspond to the emission of four bottom 
quarks, are infrared finite since the infrared singularities are softened into log- 
arithms of the bottom quark mass and can be calculated in D = 4 dimensions 
and integrated numerically without any trouble. 

The infrared divergences of the transition probabilities Cll and C12 that 
correspond to the emission of a pair of light quarks can be regulated by introduc- 
ing a small light quark mass. Nevertheless, we chose Dimensional Regularization. 
Since kinematically massless quarks and gluons are equivalent we can integrate 
the B77, B87, Cll and C12 transition probabilities in the same system 3-4. 

For the transition amplitudes B77 and B87 we have to solve the following 
basic phase space integral 

PSU)— . (4.47) 

2/34 



As we said we work in the system 3-4, (B.15). All the integrals factorize giving 

dPS(3) (2/134, 2/234) • 



rise to a simple result in the eikonal region 

1 



dP5(4)- 1 S (47r) 



2/34 167r 2 T(l - e) 



log w 



(4.48) 



As we have mentioned, the transition probabilities B77 and B87 contain in 
principle double infrared poles. Nevertheless, the sum of them, i.e., the square 
of the sum of amplitudes B7 and B8, is proportional to 

/ Vl3_ _ 2/23 \ 1 
V 2/134 2/234 / 2/34 

Expanding for small values of 2/34 we find 

A ft \ 2 a , ( 1 - 2r b - 2/134 - 2/234 T b T h \ 1 _ 1/2 , 

4«(1 - v) cos 5 — — + 0(2/ 34 ) . 

V 2/1342/234 2/134 2/ 2 34 / 2/34 (459) 



(4.49) 
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T4(a) 



T0(b) 




'f)-<L<L<b * 





T6 



T0(a) 






Bl 



B7 






B2 



B8 



Figure 1 . Cuts over a bubble diagram with three gluon vertex 
and transition probabilities with related IR structures in the 
three-jet region. The other possible cuts lead to two-jet final 
states. 



Therefore, only simple poles can appear. Same arguments apply for the transi- 
tion probabilities Cll and C12, that we grouped under Class G. 

This completes the calculation for the collincar divergent diagrams grouped 
in Class C and Class G. Intermediate steeps are rather involved but the final 



2. COLLINEAR DIVERGENCES 



49 



result can be written in a simple manner in terms of the Born amplitude 

1 



iVClass C) = ±N c as (47r)£ 



2! 4tt r(l - e) 2m.z 



w 



dPS(3)(y 13 ,y 23 ) 



e 

10 w 



— +2 + I w {y 23 ) + I w {y l3 ) 

| T0(a) +T0(b) | 2 , 



(4.51) 



3 e 



G2 



where 



I w {x) 



1 1 

2? + Te 



n , 1 , 2 ^ 1 i 2 r &W _ ^_ 

2 4 



log -1 + _ log z — - - 1ok z — - — + O ( —] 
2 9, lufe x 2 ^ 4 iu S 3,2 2 & x 2 A \ x 2 J 



(4.52) 



and 



I\„ (Class G) = T R (N F - 1 



a s (4tt) £ 



47rr(l - e) 2m z 
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with Tfl = 1/2 and (Np — 1) the number of light quarks. 

Since, as we mentioned, p 3 4 behaves in the IR region as the momentum of 
a pseudo on-shell massless particle, p| 4 — > 0, we have relabelled the integration 
variables (2/134, 2/234) as (2/13,2/23)- First contribution of Eq. ( |4.5l| ) comes from 
the transition probabilities B71, B72 and B82 and cancels the IR divergences 
of the one-loop diagrams T4 and T6. To show it explicitly needs a consider- 
able amount of algebra. We get terms proportional to the simple infrared pole 
and terms proportional to the C03 one-loop function that exactly cancels the 
I w function divergent piece. Nevertheless, we can justify such cancellation us- 
ing simple graphical arguments as we have shown in hgure [l| In fact, this was 
the main argument we used to classify all the transition probabilities in differ- 
ent groups such that IR divergences cancel independently. Second contribution 
comes from the transition amplitudes 577 and B87 and together with the Class 
G contribution is cancelled against 
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7_b + log47r, coming from the gluon wave function renormal- 



INFRARED CANCELLATIONS 



CHAPTER 5 



Numerical results and conclusions 



We have seen in the previous chapter how infrared singularities cancel be- 
tween the three parton one-loop and the four parton tree-level diagrams con- 
tributing at NLO to the Z-boson three-jet decay rate. In fact, due to the Bloch- 
Nordsiek p7[ ] and Kinoshita-Lee-Nauenberg 83 , 96] theorems, we knew from the 
beginning they should do. The cancellation of the IR divergences is of course the 
first test of our calculation. Nevertheless, the real challenge is in the calculation 
of the finite contributions. 

Ultraviolet singularities arise in loop integrals once the integration over the 
internal momentum is performed and in the form of a singular Gamma function. 
Moreover, we are allowed to expand in e before integrating over the Feynman pa- 
rameters. Main problem of IR divergences is that they appear always at the bor- 
ders. At the border of the one unit side box defined by the Feynman parameters 
in loop integrals, and at the border of the phase space for the tree-level transi- 
tion amplitudes. Therefore, we are forced to perform in arbitrary D-dimcnsions 
the full calculation in the infrared region. This makes the computation of their 
finite contributions extremely difficult. Even for the finite contributions strong 
cancellations occur among different groups of diagrams making very difficult the 
numerical approach. 

We have taken as guide line the massless result [23, 55, 91| although the 
IR structure in the massive case is completely different from the massless one. 
With massive quarks we loose all the quark-gluon collinear divergences. The 
amplitudes behave better in the IR region. The disadvantage however is the 
mass itself because we have to perform quite more complicated loop and phase 
space integrals. Furthermore, the gluon-gluon collinear divergences, leading to 
IR double poles and whose finite contributions are harder to solve, are still 
present. 

Our procedure was as follows. For the one- loop transition probabilities we 
were able to reduce the loop integrals to a few scalar one-loop n-point functions 
whose finite parts are under control, see Appendix |a|. We drop all the terms 
proportional to 1/e 2 and 1/e since, as we have seen in chapter ||, they cancel 
against the four parton real corrections to the three-jet decay rate. Since the 
boundary of the three-body phase space, including cuts, is a complicated func- 
tion we perform the remaining phase space integration over the one-loop finite 
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contributions with the help of VEGAS [^9, IOC], an adaptative Monte Carlo 



FORTRAN algorithm for multidimensional integration. 

We have seen that the four parton transition amplitudes have a simple and 
well known behaviour in the soft and collinear regions, basically proportional 
to the lowest order three parton transition probabilities, making the integration 
over these regions of the four-body phase space feasible analytically. Our strat- 
egy was to exclude from the numerical integration domain the singular regions. 
We impose a cut, we called w, with w very small, in the energy of the soft glu- 
ons. For the diagrams with collinear singularities we cut on the two momenta 
invariant variable t/34 = 2(p 3 • Pi)/s. In the excluded phase space regions we 
replace the transitions amplitudes by their limiting values and perform analyti- 
cally the integration over this phase space singular region in arbitrary D = 4 — 2e 
dimensions (we only need to analytically integrate three of the five independent 
variables defining the four-body phase space). The analytic integration gives 
the IR poles in 1/e 2 and 1/e that again we drop[] and some finite contributions 
that we still have to integrate numerically over an effective three-body phase 
space. Finally, the full transition probabilities, calculated in D = 4 dimensions, 
are integrated numerically in the three-jet region of the four-body phase space 
above the w cut. 

Both results, the semi-analytical one and the full-numerical one depend on 
the cut w. First requirement should be of course the sum of them is independent 
of w. We have to keep in mind this is an approximate method since we have sub- 
stituted in the singular region the transition amplitudes by their limiting values. 
The smaller value of w the better approximation we get. On the other hand, 
since the integrand near the boundaries of the singular region increases steeply, 
typically like logiu, the w cut can not be very small otherwise we immediately 
loose accuracy in the numerical integration. In practice we should find a compro- 
mise over these two conflicting requirements to achieve the best efficiency. We 
have studied the w-dependence of the final result and we found w/mz ~ 10 -5 
provides a good approximation. 

This method is called in the literature |90| the phase space slicing method 
as we exclude from the numerical integration domain a slice of the phase space. 
It was applied for instance by [11, 12,|66, 57, 72|. Another method of analytic 



cancellati on of I R si ngu larit ies, calle d the subtraction method, has been devel- 



oped by |34 , 55] , p6| , |58| , pi] , 92 , |101| . The idea is to subtract and add back a 
quantity that must be a proper approximation to the differential real cross sec- 
tion such that it has the same pointwise singular behaviour (in D-dimensions) 



1 The IR poles appear multiplied by the Born amplitude calculated in D-dimensions and by 
some e-dependent factors like (47r) e /r(l— e) that in principle would generate finite contributions 
for e — » 0. Nevertheless, since both the one-loop and the four parton transition probabilities 
share the same factors they will cancel as well as the IR poles and therefore we are allowed to 
drop these extra contributions. 
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and is analytically integrable (in D-dimensions) over the one parton subspace 
leading to the soft and collinear divergences. 



1. Tests 



As we mentioned the introduction of massive quarks complicate extremely 
the calculation of the three-jet decay rate in spite of the fact that some of the 
singularities are softened into logarithms of the quark mass. Therefore, to check 
our results against already tested simpler calculations is a crucial point. 

The knowledge of the transition amplitudes is a well established result in 



the massless case [55, 57, 87] as well as the three- and four-jet fractions in the 
JADE, E and DURHAM jet clustering algorithms ||23L^l|. The four-jet decay 



rate has been calculated for massive quarks for instance in [13, 14, 69 



First test was to check our four parton final state transition amplitudes in the 
four-jet region. Since the four-jet phase space region is free of IR singularities in 
both the massless and the massive case we are allowed to neglect the quark mass. 
The comparison with the massless result 



23 



55, 



91 



was successful as well as 
with the massive calculations [|13j, |14L |69j . Nevertheless, this is not a conclusive 
proof. First, this checks only the four parton final state transition amplitudes. 
Second, the calculation can be performed in four dimensions. And finally, since 
the four-jet region is completely IR finite, it does not matter the way we write the 
transition amplitudes, the Monte Carlo phase space integration converges quite 
quickly. The main problem is not how to calculate the transition probabilities 
but the way we write them. The ERT citeEllis81 transition amplitudes can not 
be directly used when we try to calculate the three-jet decay rate. 



1.1. Massless three-jet decay rate at next-to-leading order. Next 
steep was to test our procedure of analytic cancellation of IR divergences over 



the ERT [55 1 transition amplitudes in the three-jet region and to recalculate the 
A 1 * 1 ^ (y c ) function that defines the NLO massless correction to the three-jet decay 
rate. Furthermore, this test allowed us to check our phase space, algorithm cuts 
and numerical integration routines in the massless limit. We remind the three-jet 
decay rate can be written as 



r 3i = m z^TT~A S -i9v H v(yc,n) + g 2 A H A (y c ,r b )}, (5.1) 



where r b = rn^/in%, y c is the resolution parameter, and gv and g A are the vector 
and the axial- vector neutral current quark couplings. At tree-level and for the 
bottom quark gv = — 1 + 4s^/3 and g A — 1 in the Standard Model. In the 
massless limit we wrote 



H viA) =A^(y c ) + ^A^(y c ), (5.2) 
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Table 1. Results of the five parameter fits to the function 
A (1) (y c ) = En=o fc ™ l °S n Vc in the range 0.01 < y c < 0.20 



Algorithm 



k 



ki 



A: 2 



EM 17.25(12) 20.37(10) -7.320(47) 

JADE 17.872(26) 26.599(21) 1.265(10) 

E 12.41(13) 20.22(10) -1.912(49) 

DURHAM -.021(13) Jf ^Z fl ffl)) ft ftfffif^ 



fc 3 

-12.621(21) 
-9.5644(48) 
-11.173(23) 
-3.3720(25) 



-1.9200(45) 
-1.8547(10) 
-1.7920(47) 
-.46956(49) 





Figure 1. The function which determines the massless 
NLO three-jet decay rate for the EM, JADE, E and DURHAM 
jet clustering algorithms. 



1. TESTS 



55 



where we have taken into account that for massless quarks vector and axial 
contributions are identical] 



The EM jet clustering algorithm was introduced for the first time in 25 1 be- 
cause in this algorithm at LO the two- and three-jet decay rates can be calculated 
analytically for massive quarks. The function A^> (y c ) is therefore unknown in 
this case. We have calculated it for the first time. Furthermore, we have not only 
tested but tried to improve the accuracy in the knowledge of the A^(y c ) func- 
tion for the JADE and DURHAM algorithms. Quark mass corrections are very 
small in these two jet clustering algorithms, specially in DURHAM, and a better 
determination of the A^(y c ) function contributes to a better understanding of 
the quark mass effects. 



To simplify the use of our results and to compare them with the [23, 91] 
result we have performed a simple five parameter logy c power series fit to the 
A (Uc) function for the different jet clustering algorithms 



4 

A^(y c ) = Yk n \og n y c . (5.3) 



The results of our fits for the range 0.01 < y c < 0.20 are presented in table We 
plot in figure ^ the A^ (y c ) function for the four algorithms we have considered. 

1.2. The massless limit test. Finally, we integrate in the three-jet region 
the massive quark transition amplitudes we have calculated and we try to see if 
extrapolating our result to the massless limit we can reach the ^l*- 1 ^ (y c ) function 
of the previous section. This is the main test of our calculation. 

In figures || and [| we present our result for the vectorial contribution to 
the 0(a 2 s ) three-jet decay rate of the Z-boson into bottom quarks for the JADE 
and E algorithms. We have performed the calculation for different values of the 
bottom quark mass from 1 to 5(GeV) for fixed y c . We remind we factorized the 
leading dependence on the quark mass in the functions -Hy(A) &s follows 

H V(A) = AWfofe) + r b B§\ A) {y c ,r b ) + ^ (a^( Vc ) + r b B$ A) (y c ,r b )) . 

(5.4) 

In the JADE algorithm we can see that for big values of y c the NLO correc- 
tions due to the quark mass are very small and below the massless result. Notice 
they increase quite a lot for small values of y c and give a positive correction that 
will produce a change in the slope of the LO prediction for In any case 

we recover the massless limit and a linear parametrization in the quark mass 
squared would provide a good description. The DURHAM and EM algorithm 
would exhibit a very similar behaviour to the JADE algorithm. 

The E algorithm behaves also linearly in the quark mass squared although 
only for big values of y c . Corrections in the E algorithm are always very strong. 



2 We do not conside r the small 0(oP s ) triangle anomaly jr3|. A^(y c ) = B(y c )/4 with 
B{y c ) defined in §31 0. 
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FIGURE 2. NLO contribution to the three-jet decay rate of Z —> 
bb for bottom quark masses from 1 to h{GeV) and fixed y c in 
the Jade algorithm. Big circle is the massless case. 
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FIGURE 3. NLO contribution to the three-jet decay rate of Z — > 
bb for bottom quark masses from 1 to 5(GeV) and fixed y c in 
the E algorithm. Big circle is the massless case. 



The reason is the following, the resolution parameter for the E algorithm ex- 
plicitly incorporates the quark mass, yij = (pi +pj)/s, i.e., for the same value 
of y c we are closer to the two-jet IR region and the difference from the other 
algorithms is precisely the quark mass. This phenomenon is already manifested 
at the LO. For massive quarks the behaviour of the E algorithm is completely 
different from the other algorithms, ft is difficult to believe in the E algorithm 
as a good prescription for physical applications since mass corrections as so big. 
However for the same reason, it seems to be the best one for testing massive 
calculations. 
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Furthermore, since we have encoded the ERT formulas, we were able to 
perform some partial tests, not only the final three-jet decay rate. The IR 
structure of the massless transition amplitudes of [55] is completely different 
from the massive case. This means, we can not test our calculation amplitude 
by amplitude. Nevertheless, in the previous chapters we were able to split our 
transition probabilities into different groups such that IR singularities cancel 
independently. Checking can be performed group by group. The difference lies 
in the singularities that are softened into logr;, for massive quarks. They must 
cancel too. In the massless case we have to cancel these singularities analytically. 
In the massive case we do it numerically. We have not calculated analytically 
the terms proportional to logr^. A future improvement of our calculation could 
be to do that. 



2. Massive three-jet decay rate at next-to-leading order 

To conclude, we present our final calculation for the By} A Jy c ,rb) functions 

defining the NLO mass correction to the three-jet ratio R^ d in the four jet clus- 
tering algorithms EM, JADE, E and DURHAM, see figures f§ |, | and g. We 
have performed the calculation for mi, — 3(GeU) and mj, = b(GeV) pole masses. 
Error bars come from the Monte Carlo phase space numerical integration. As be- 
fore, to simplify the use of our results, we have performed simple five parameter 
logy c power series fits to the ratio of the function By A (y c ,rb) and the massless 
NLO function A^\y c ). We use the following parametrization: 

4 

B^ A (y c ,r b )/A^(y c ) = £ k VA log n y c . (5.5) 

The results of our fits for the range 0.01 < y c < 0.20 (0.025 < y c < 0.20 for 
JADE and E algorithms) are presented in table ^|. For the moment we perform 
independent fits to the mj = 3(GeV) and mb = h(GeV) results. 

DURHAM seems to be the most promising jet algorithm to keep under 
control the higher order corrections since not only the NLO massless function 
A^\y c ) is the smallest but also because the ratio By\(y Cl rb)/A^ (y c ) holding 

the mass information is the smallest too. We get By A (y c ,r b )/A^(y c ) ~ —10. 
On the other hand, mass corrections are so big in the E algorithm, for the 
reasons we mentioned in the previous section, and LO and NLO predictions for 
R^ d are so different that perturbation theory probably breaks down in this case. 
The normalized mass correction in the JADE and EM algorithms are around 
two and three times the mass correction in the DURHAM algorithm, we obtain 
By\(y Cl ri,)/A^ (y c ) ~ —20,-30 respectively, although still under reasonable 
values. 

In all the three algorithms, EM, JADE and DURHAM, the B$ A (y c , r b ) /A& (y c ) 
ratio is negative since massive quarks are expected to radiate less gluons than 
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Table 2. Results of the five parameter fits to the function 
B^ A (y c )/AW(y c ) = £* =0 k^ A \og n y c in the range 0.01 < y c < 
0.20 for EM and DURHAM algorithms and 0.025 < y c < 0.20 
for JADE and E for bottom quark pole masses of 3(GeV) and 
5(GeV). 



Algorithm 




fci 


fc 2 


fc 3 


&4 


EM 


k v (3GeV) 


6.9(48) 


90.2(30) 


84.15(93) 


27.79(27) 


2.957(47) 


k A 


-54.5(16) 


-43.32(94) 


-6.85(28) 


3.024(84) 


0.664(14) 


k v (5GeV) 


-10.7(52) 


65.4(31) 


71.98(93) 


25.39(27) 


2.798(47) 


k A 


-38.4(16) 


-24.09(95) 


1.49(28) 


4.731(83) 


0.800(14) 


JADE 













k v (3GeV) 396.3(39) 676.8(25) 410.59(92) 109.93(32) 10.941(75) 

k A 85.5(12) 170.40(75) 111.77(28) 33.256(97) 3.747(23) 

k v (5GeV) 224.2(38) 403.2(25) 249.43(92) 68.35(32) 6.989(75) 

k A 94.0(12) 181.87(75) 117.55(28) 34.689(97) 3.893(23) 



E 



k v (3GeV) 


589.6(52) 


1080.1(32) 


759.5(12) 


234.24(40) 


27.599(92) 


k A 


-584.7(17) 


-1059.1(10) 


-673.70(38) 


-183.63(13) 


-17.305(30) 


k v (5GeV) 


621.8(52) 


1112.7(32) 


767.8(12) 


233.85(40) 


27.343(92) 


k A 


-535.4(17) 


-975.3(10) 


-622.80(38) 


-170.11(13) 


-16.004(30) 


DURHAM 












k v (3GeV) 


105.0(38) 


183.4(25) 


107.27(77) 


27.17(23) 


2.498(38) 


k A 


69.8(11) 


110.34(76) 


61.53(23) 


15.364(69) 


1.421(11) 


k v (5GeV) 


132.5(38) 


212.5(25) 


120.40(77) 


29.99(23) 


2.739(38) 


k A 


74.9(11) 


116.88(76) 


65.51(23) 


16.435(69) 


1.525(11) 
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massless quarks. Nevertheless, for enough small value of the resolution param- 
eter y c , depending on the algorithm, the mass correction becomes positive and 
steeply increasing. This means there is one point where the NLO mass correc- 
tions are zero. 

This phenomenon happens because for enough small values of y c the four 
parton transition probabilities contribute mainly to the four-jet fraction and 
therefore in the three-jet region the virtual corrections dominate. This translates 
into a change in the slope of the LO prediction for R^ d , i.e., the NLO prediction 
for R^ d going down for small values of y c becomes increasing and closer to unity. 
For the range of the resolution parameter y c we are considering, 0.01 < y c < 0.20, 
we can see in figures [| and |^ this already occurs in the EM and JADE algorithm 
for y c ~ 0.02. For the DURHAM algorithm we have to wait for a resolution 
parameter smaller than y c < 0.01. In fact, although this phenomenon is a sign 
that perturbation theory breaks down and therefore our perturbative prediction 
is meaningless, we can see in figure ^ that for the preliminary DELPHI results 
in the DURHAM algorithm this tendency really happens. 



2.1. Estimate of higher order contributions. Up to now we made all 
the calculations with the assumption /i 2 = s = m 2 z . This does not restrict 
the validity of our result since the /i dependence is completely determined by 
the Renormalization Group. The scale dependence is introduced in R^ d by the 
replacements 



n 



fb(p) { 1 + 2 



,a s (/x) 



4 . ml(fi) 



3~ l0g - fj, 2 
>(o) 



(5.6) 



Assuming we have factorized in the functions By' A (y c ,ri,) and By' A (y c ,rb) the 
leading dependence on the quark mass and the remnant mass dependence is 
soft 0, i.e., By^ A (y c ,r b ) ~ By^ A {y c ) the three-jet fraction ratio reads 



R 




log- 



A(°Hy c ) 



c v 2 



(5.7) 



3 For the moment we will assume the By A (y c ) function comes from the calculation with 
m b = 3(GeV). 
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At the lowest order we define 



R\ d 



cv 



ggV) 



CA 



(5.8) 



with m° = mzyrj. We can solve iteratively Eq. (5.7) to get 



n(m z ) = r\ 



47o 

, a s (m z )\ /3 



a 



Rf - 1 



cv 



+CA 



-YlCA 



Bf{y c ) 
4 



X 0) 0/ c ) 

gj 1} (yc) 



log- 



21og^ 



--21og^ 
.3 M 



A(°)(y c ) 



(5.9) 



12c V 



where /3o = 11 — 2/3Np and 70 = 2, with the number of flavours Np = 5. 

Using our fit for a pole mass of mj, = 3(GeU) and following Eq. (5.7) we have 
plotted in figures ^, ||; |^ and ^ our prediction for R^ d for a running bottom quark 
mass at the mz scale of rhb(mz) = 3(GeV). We include also our prediction for 
/i = mz/2 and fj, = 2 for fhf,(mz) = 3(GeV) fixed that could be taken as 
our estimate of the higher order uncertainty. 

Conversely, for a fixed value of the resolution parameter y c and a fixed value 
of the ratio R^ d we have explored the /i dependence of our prediction for the 



bottom quark mass at the Z-boson mass scale, fhb(mz) , following Eq. (5.9), see 
figures [|, ||, H and 0. We include, to compare, the LO /i-dependence. 

Taking as the theoretical uncertainty for the running bottom quark mass 
at mz half of the difference from the = mz/2 and (i = 2 mz predictions 
we find as expected the smallest uncertainty correspond to the DURHAM algo- 
rithm, Ami, — HO(MeV). JADE and EM share the same behaviour, Am& ~ 
130(MeU), slightly bigger than in the DURHAM case. Finally, the E algorithm 
grows to Amj — 220(MeV). Case we enlarge the lower limit on /j, we get for 
fi = m z / 10, Am b ~ 250, 300 and 500(MeU) respectively. 



3. Discussion and conclusions 

In this thesis we have analyzed at NLO the quark mass dependence of the 
three-jet observable 

gM _ £W/E? 

3 rUy c )/r d ' 
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The EM algorithm at NLO 




Figure 4. EM algorithm: functions By{y c ) and B\'(y c ) nor- 
malized to he massless NLO function A^(y c ) in the range 
0.01 < y c < 0.20, three-jet fraction ratio R^ d (solid lines arc 
LO for m b = 3(GeV) and m b — 5(GeV), short dashed line is 
the NLO prediction for a running mass fhb(mz) = 3(GeV) and 
long dashed lines for /j, = mz/2 and fi — 2 mz) and NLO scale 
dependence for the bottom running mass extracted from a fixed 
value of the ratio R^ d (dashed line is LO). 



for different jet clustering algorithms such as EM, JADE, E and DURHAM. In 
particular, we have studied the possibility of extracting the bottom quark mass 
from R\ d at LEP. 
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The JADE algorithm at NLO 




Figure 5. JADE algorithm: functions By\y c ) and B^(y e ) 
normalized to he massless NLO function A^(y c ) in the range 
0.01 < y c < 0.20, three-jet fraction ratio R^ d (solid lines are 
LO for nib — 3(GeV) and mt, = 5(GeV), short dashed line is 
the NLO prediction for a running mass fhf,(mz) = 3(GeV) and 
long dashed lines for \i = mz/2 and ju = 2 mz) and NLO scale 
dependence for the bottom running mass extracted from a fixed 
value of the ratio R\ d (dashed line is LO). 



Since quarks are not free particles the study of its mass needs a precise 
theoretical framework and, in fact, they should be treated more like coupling 
constants. The two most commonly used quark mass definitions are the pertur- 
bative pole mass and the MS scheme running mass. We reviewed in Chapter |l| 
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The E algorithm at NLO 
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Figure 6. E algorithm: functions By (y c ) and B { A '(y c ) nor- 
malized to he massless NLO function A^(y c ) in the range 
0.01 < y c < 0.20, three-jet fraction ratio R^ d (solid lines arc 
LO for m b = 3(GeV) and m b — 5(GeV), short dashed line is 
the NLO prediction for a running mass fhb(mz) = 3(GeV) and 
long dashed lines for \i = mz/2 and fi — 2 mz) and NLO scale 
dependence for the bottom running mass extracted from a fixed 
value of the ratio i?g d (dashed line is LO). 



some of the most recent determinations of all the quark masses, in particular, the 
bottom quark mass extracted from the bottomia spectrum from QCD Sum Rules 
and Lattice calculations. The bottom quark perturbative pole mass appeared 
to be around M& = 4.6 — 4.7(GeV r ) whereas the running mass at the running 



64 



5. NUMERICAL RESULTS AND CONCLUSIONS 



The DURHAM algorithm at NLO 




Figure 7. DURHAM algorithm: functions By (y c ) and 
B^ {He) normalized to he massless NLO function (y c ) in the 
range 0.01 < y c < 0.20, three-jet fraction ratio R^ d (solid lines 
are LO for nib — 3(GeV) and mt = 5(GeV), short dashed line 
is the NLO prediction for a running mass fhb(m,z) = 3(GeV) 
and long dashed lines for \i = m z /2 and fi = 2 m z ) and NLO 
scale dependence for the bottom running mass extracted from 
a fixed value of the ratio i?g d (dashed line is LO). 



mass scale read ffib{rhb) = (4.33 ± 0.06)GeV. We analized its running until the 
Z-boson mass scale and we found fhb(mz) — (3.00 ± 0.12)GeV. 

Since for the bottom quark the difference between the perturbative pole mass 
and the running mass at the mz scale is quite significant it is crucial to specify, 
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Figure 8. Preliminary experimental results for the ratio R^ d 
from DELPHI Collaboration. 



in any theoretical perturbative prediction at mz, which mass should we use. In 
practice, this means we have to go up to NLO, where a renormalization scheme 
has to be chosen, if we want to extract a meaningful value for the bottom quark 
mass, since the LO does not allow to distinguish among the different quark mass 
definitions. 

We have explored a new method: the study of three-jet observables at LEP, 
different from QCD Sum Rules and Lattice calculations, for determining the 
bottom quark mass. Among other advantages, it would allow to extract the 
bottom quark mass far from threshold, in contrast to the other methods described 
above. 

From our calculation, we have estimated, by exploring the /i-dependence 
of our prediction at mz, a theoretical bottom quark mass uncertainty of at 
most 250(MeU) in the DURHAM algorithm, 300(MeU) for JADE and EM and 
500(MeV) for the E algorithm. The preliminary experimental bottom quark 
mass uncertainty from the DELPHI Collaboration |6C)| , see figure ||, in the 
DURHAM algorithm is expected to be 300(MeU). Nevertheless, that is not 
the last word but the first step. The experimental error would be improved, spe- 
cially once the four experiments at LEP are taken into account, and we are still 
discussing how to improve the theoretical error. Resummed expressions would 
probably give a better result. 
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Therefore, even if, for the moment, this method can not be considered com- 
petitive with QCD Sum Rules and Lattice predictions, a total error of 300(MeV) 
in the DURHAM algorithm, slightly bigger in JADE and EM algorithm, should 
be enough to check, for the first time, how the bottom quark mass evolves, from 
low scales /i ~ h(GeV) to the mz scale, following the Renormalization Group 
prediction, in the same way the running of the strong coupling constant has been 
checked before. 

Among other applications, our calculation would provide a better under- 
standing of the systematic errors in the measurement of a s (mz) from 66-production 
at LEP from the ratio of three to two jets [ 130| ] and therefore a better measure- 
ment of the strong coupling constant a b s (mz). 

Furthermore, quark mass effects are expected to be quite more important in 
top production at NLC (Next Linear Collider) and our calculation would be of 
a great interest. 

In this thesis, we have concentrated in the study of the three-jet observable 
i?3 d ratio. To do that, we had to calculate the NLO matrix elements for the 
decay of the Z-boson into massive quarks. Main problem of such calculation 
was, in addition to UV divergences, the appearance of IR singularities. We were 
able first, by simple graphical arguments, to classify our transition amplitudes 
into different groups such that IR singularities cancel independently and second 
to build matrix elements free of singularities since we exclude from the phase 
space the IR region where integration is performed analytically. Once the IR 
behaviour is under control, it would be possible to build, with our matrix ele- 
ments, a Monte Carlo generator opening the door to the study of other kind of 
three-jet observables (angular distributions, differential cross sections, etc . . . ) 
where massive quarks are involved. Furthermore, the two-jet decay rate can be 
extracted by subtracting from the well known total decay rate the three- and 
the, full infrared finite, four-jet decay rates. Therefore, the study of two-jet 
observables is also possible from our calculation. 

Concluding, in a previous paper [ 25 ] , we raised the question of the possibility 
of measuring the 6-quark mass at LEP by using three-jet observables. In close 
collaboration with the DELPHI experimental group at Valencia J60| |, we have 
developed the necessary theoretical tools in order to extract a meaningful value 
for the 6-quark mass from LEP data, in particular, the NLO calculation of the 
three-jet decay rate of the Z-boson into massive quarks. Clearly, more work 
has to be done but the effort is worth since it will allow for an independent 
measurement of nib, far from threshold, and it will provide a test of the QCD 
Renormalization Group predictions. 



APPENDIX A 



Loop integrals 



We present in this appendix all the scalar one loop integrals we encountered 
in the calculation of the diagrams of figure |l|. Since we are just interested in the 
interference transition probabilities of these diagrams with the born amplitudes 
we will only need the real part of such n-point scalar one loop integrals. Although 
we don't specify it, the following results refer, in almost all the cases, just to the 
real part of these functions. 



They are the source of UV divergences. 

For the scalar one-point integral in D = 4 — 2e dimensions we have the 
following result 



1. One- and two-point functions 



One- and two-point one loop functions are defined in D dimensions as 




(A.l) 



A(m) = m 2 (A - log — 
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1), 



(A.2) 



with 



A = 7 + log 4-7T . 



(A.3) 
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the UV pole while for the two point functions we get 
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(A.4) 



where 



r b = m 2 b /m 2 z , y^ = 2{ Pi ■ Pj)/m 2 z 



P = VT Z ^ r b, 012 = Jl- 4 Jo" > ( A ' 5 ) 

2/12 + 2r 6 



1-/3 1-/3; 



Cl2 



12 



1 + /? 1 + ^12 

We have set to zero the i? (0,0,0) integral. Strictly speaking that is not 
completely true. On analyzing this integral in detail we encounter it receives 
contributions in terms of UV and IR divergences separately, i.e., we should write 

flo(0,0,0) = A D / 7^7^ = — - — ■ (A.6) 
J (27r) u A: 4 e uv e IR 

If we try to distinguish between UV and IR poles we have to keep them. Since 
this kind of integrals just change the label of the divergences but not the number 
of poles in 1/e we will identify em — euv- We can try to follow the origin of 
the UV and IR poles but we don't gain nothing. Ultraviolet singularities are 
cancelled by the renormalization procedure that is very well known. 

2. Three-point functions 

The scalar three-point integral reads 

Co(A, B, C, m\, ml, m 2 ) = 



16tt 2 



^ ° I (2tt) d [k 2 - m 2 } [{k + P i) 2 - ml] \{k + p 2 ) 2 - m 2 } (A ' ?) 



2. THREE-POINT FUNCTIONS 



(i!) 



where A = pf, B — (pi — P2) 2 and C = p\. After Feynman parametrization and 
momentum integration we get 



C {A,B,C,ml,mlml) 



r 3 
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dx / dy 



(ax 2 + by 2 + cxy + dx + ey + /) 3 



(A.; 



where 



a = C , 
b = B , 

c = A — B — C , 
d = m 2 — ni\ — C 

4 

f = m\ — irj . 



(A.9) 



e = m\ — m% + C — A 



For A, B and C different from zero and if the scalar three-point function 
does not contain infrared divergences the integral can be performed in D = 4 



dimensions and can be expressed as a sum over twelve dilogarithms ]77 

3 2 



Xi yij 
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with the Spence function or dilogarithm function defined as [ 113( | 

log(l-f) 



dt- 
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d + 2a+ (c + e)a 
d + ea 



(A.10) 



(A.ll) 



(l-a)VA ' 
d + ea 



-(c + e) ± y/(c + ef - 4b(a + d + f) 



(A.12) 



2b 



-(d + e) ± y/{d + e) 2 - 4/(a + b + c) 
2(a + 6 + c) 

-d± v^ 2 - 4a / 
2a ' 
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where a is a real solution of the equation 

ba 2 + ca + a = 

and y/\ = c + 2ab. 

We define the following one loop integrals 



I6n' 
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(A.13) 



167T 



' C02(y 13 )= M 4 - D ' '' k 



(2tt) d k 2 [(k + Pl3 ) 2 - m 2 ][(fc + Pl f - ml 



r C03(yi 3 ) = 



167P 



4-Z5 



:C04(j/i 2 ) = /x 



4-D 



1 



a 7 



(2tt) d fc 2 (fc - p 3 ) 2 [(k + Vl f - m; 

d D /c 1 
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167T 



4-D 



d D k 



(2ir) D k 2 [(k + Pl f - to 2 ] [(jfe - p 2 ) 2 - m 2 b ] ' 



For integral C01 we can apply directly Eq. (A. 10) and take just its real 
part. For C02 and C04 we can calculate them either directly from the Feynman 
parametrization (A. 8) or from Eq. (A.1C) by taking the appropriate limit. At 
the end we get a quite simplified result [J 
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1 - {y 12 + 2r b ) 2 

with A = rh/(rft + yi 3 ). The three-point functions C03 and C05 a re IR divergent 
and therefore we can not apply the general result of Eq. ( A.1C ). In these two 
cases we obtain, after integration in D = 4 — 2e dimensions n 
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(A.16) 



*It is understood only the real part. 

2 C05 is the same one-loop integral that appears in the lowest order one-loop diagram. 
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C05(y 12 ) 



(4tt) 



T(l - e) (j/12 + 2r b )(3 l2 



logr b logci2 - 2L(y 12 ) 



where 



7T 2 1 

L(yi2) = Li 2 (c 12 ) + — + log(l - ci 2 ) log C12 - 1 log 2 C12 . 



(A.17) 



(A.18) 



3. Four-point functions 

In our one loop matrix elements we encounter the following four-point func- 
tions 



4-D 



// 



d D k 



(2Tt)Dk*[{k + Pl y-mi][{k + p l3 y-mi][{k-p 2 y- m i] ' ! U!,J 



167T 2 



£>06(yi 3 , 2/23) 
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(27r)^A : 2 (fc+p 3 ) 2 [(fc+p 13 ) 2 -m 2 ][(fc- i , i - ,„,, 



p 2 ) 2 - TO 21 



(A.20) 



Both four-point functions are IR divergent. The integral D05 has simple 
poles in 1/e while D06 presents double poles. To get their finite part is not an 
easy task but in order to extract the divergent piece we just need to split its non 
on-shell propagator. 

For D05 we can write 

1 _ 1 k 2 + 2k ■ p 13 

(fc+pi 3 ) 2 -m 2 Ph-m 2 b (p 2 13 - m 2 )[{k + p 13 ) 2 - m 2 } ' ^ ^ 

The IR divergence is isolated in the first term of the righthand side of the previous 
equation and give rise to C05(yi 2 )/yi3. The other term generates two IR finite 
one loop integrals. 

For D06 things are more complex. The splitting of the non on-shell propa- 
gator does not give directly the divergent piece due to the presence of IR double 
poles. First of all we have to notice D06 is invariant under the interchange of 
particles 1 and 2, i.e., the result of the integral should be symmetric in the two 
momenta invariants j/13 and ^23- We found the proper way for extracting the IR 
piece is as follows 

D06{y 13 ,y 23 ) = — C03(y 23 ) + — C03(y 13 ) + finite terms . 

Via 2/23 (A22) 

We leave for a future publication the detailed calculation of the finite con- 
tributions of the box one-loop integrals D05 and D06. 
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4. Explicit calculation of C03 



For the C03 function a Feynman parametrization, (A. 8), with both Feynman 
parameters running from to f is more convenient because the two integrals 
decouple each other 

C03 = -(4tt)T(1 + e) f aT 1 " 2 ^ 
Jo 

[(n + Vi3)y 2 - (2r fc + y 13 )y + r b ~ irj] 1 £ dy . (A.23) 



The integral over the x-parameter is immediate and we left just with 



C03 = (4^r r(1 2 ^~ £) (n + yu)- 1 -' 



{l-y)- l - f {yv-y)-^dy 



(A.24) 



where yo is the root of the y-parameter polynomial. The last integral gives rise 
to an Hypergeometric function 
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J/o 



that with the help of some mathem atica l properties ||, |ll5| can be written in 
terms of a dilogarithm function [46, 113] 
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(A.26) 



In the region of physical interest yo = r(,/(r*{, + j/13) is smaller than one and we 
should perform an analytic continuation taking into account the —irj prescription. 
The final result is quoted in Eq. (A. 16). 



5. Passarino-Veltman reduction 

We illustrate here a simple example of how to reduce vectorial and tensorial 
loop integrals |l!4| to the scalar one loop functions defined in the previous 
sections. Let's consider the following set of integrals 

i—D f dPk fc^; k^k u i 

^ J {2^k 2 [(k+ Pl Y~ml][{k-p 2 f-ml] ~ 16^° W 
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By Lorentz invariance 

C0» v = g^C 20 + tftfCii + P^P 2 C 22 



(A.28) 



Now, we have to perform all the possible contractions with external momenta to 
extract the coefficients Ca . 



Let's try with 
d D k 



,4-D 



(k-pi) 



(2rr)^ fc 2 [(fc + Pl y - ml] [{k - P2 ) 2 - m 2 b ] 



167T 1 



[pfCn + (px ■ p 2 )C 12 ] . (A.29) 



The scalar product in the numerator can be written, including the on-shell con- 
dition p\ = ml, as 



2(k ■ Pl ) = (k+ Pl ) 2 -k 2 -pj = [(fc + Pl f - ml] - k 2 



(A.30) 



Therefore, our initial vectorial three propagators integral reduces automatically 
to the sum of two scalar two-point functions 



,4-D 
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167T 
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(A.32) 



Same argument can be applied to the scalar product with P2 . In this simple 
case the final result is as follows 



C n 



C\ 2 — 



1 



2/12 - 2r & 



[B (p 2 i 2 ,ml,ml) - J3 (m 2 ,0,m 2 )] 



(A.33) 



To reduce the tensorial integrals to scalar integrals we have to follow the same 
mathematical procedure. Solutions appear to be quite more complex, specially 
when we have to treat with box integrals, i.e., one loop integrals with four 
propagators. To solve all the P asarino-Veltman reductions we needed we have 
used the FeynCalc 1.0 [45, 107| algebraic package for Mathematica 2.0. whereas 
all the Dirac algebra was performed with HIP 
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6. The problem of 75 in D-dimensions 

Multi-loop calculations with dimensional regularization often encounter the 
question of how to treat 75 in D dimensions. Occasionally the problem can 
be circumvented by exploiting chiral symmetry. In our case, in the limit of 
massless quarks axial and vector contributions to the three jets decay rate of the 
Z boson should be equal^J. In general, however, a consistent definition must be 
formulated. 

A rigourous choice based on the original definition of 't Hooft and Velt- 
man 



|76| 



75 = j^up^Yr ■ (A.34) 

can be found in [29|. As a consequence of the lost of anticommuntativity of 75 in 
this new definition, standard properties of the axial current as well as the Ward 
identities are violated. In particular, extra finite renormalization constant should 
be introduced to restore the correctly renormalized non-singlet axial current.. 

Nevertheless, it has been checked |30[ that for diagrams with an even number 
of 75 connected to the external current the treatment based on a naive anticom- 
mutating 75 leads to the same answer. This is our case because for unpolarized 
final state quarks only terms proportional to the square of the vector and axial- 
vector neutral current couplings survive. Therefore, in our approach we have 
restricted to work with anticommutating 75 avoiding extra complications. 



3 Singlct contribution not included. 



APPENDIX B 



Phase space in D = 4 — 2e dimensions 



The phase space for n-particles in the final state in D-dimensions [64, 102, 



103 1 (D = 4 — 2e) has the following general form 

i— l,n \ i—l,n 

- II j^-Ap* - m?)e(i3.)f D («- £ ^ (B ' 1} 

i—l,n \ i— l,n , 

Then doing several trivial integrations we have the following phase-space factor 
for the process Z — > 66 

where /3 = \/l — 4?*b with r?, = m 2 /m z , 

For the case of the decay into three particles, Z — > 6(^1)6(^2)3(^3), we get 



(B.3) 



where the function h p which gives the phase-space boundary in terms of variables 
2/13 = 2(pi ■pz)jm 1 z and 2/23 = 2(p2 -p^/rriz has the form 

fop = 2/132/23(1 - 2/13 ^ 2/23) - n>(?/i3 + 2to) 2 • (B.4) 
1. System 1-3 

For those four parton transitions amplitudes of the process Z — > 6(pi) + 
b(p2) +3(^3) +3(^*4) containing the denominator 2/13 = 2(pi ■p3)/m%, i.e., with 
gluon labelled as 3 soft, it is convenient to write the four-body phase space as a 
quasi three body decay 

q -^13 +P2+PA 

(B.5) 

^ Pi + P3 ■ 
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Figure 1. Threemomenta in the cm. frame of particles 1 and 3. 



We will refer to this as the "1-3 system" 
1 and 3 the four-momenta can be written as 
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In the cm. frame of particles 



Pi — (Ei, ... , Pi sin 9 cos 0', Pi cost 
p 2 = (E 2 , ... ,0,p 2 ) , 
P3 = Es(l, . . . , — sin#cos(9', — cost 
P4 = Ei(l, ... , sin/3, cos j3) , 



(B.6) 



where the dots inpi andj>3 indicate D — 3 unspecified, equal and opposite angles 
(in D dimensions) and D — 3 zeros in p 2 and pi. In terms of variables 



2/13 



2(pi • Ps) 
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2(P2 • P13) 



2/134 



2(P4 -Pl3) 



(B.7) 



with P13 — pi + P3, energies and threemomenta read (s = m| = 1) 



Si = 
E 2 = 
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27^ 
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2V2/13 + H. 
2/134 

2^2/13 + 



^23-4^(2/13+^) , (B.8) 
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Setting v = |(1 — cos 9), we obtain for the D-dimensional phase space in 
this system 

re ' 4 »= 16 ;w-2.) (^r' y*-*» 



*-.>^^/^ W _j£_ ( B,, 



/ -«))- £ — / ^'sin- 2e 0' 

Jo N e> Jo 



where S = 1/2! is the statistical factor, Ng> is a normalization factor determined 
such that 

[ V dO 1 sin" 26 9' = N , = 2 2e 7r r ^~ 2£) , (B.10) 
Jo r (1 - e) 

and the function 

hp = 2/1232/134(1 - 2/123 - 2/134) 

+ r b (-1 + 2(2/123 + 2/134) - (2/I23 + 42/1232/134 + 22/134)) 

+ 4r^(l -2/123 -2/134) -4rg 

+ (-l + 6r 6 -8r 2 (B.ll) 

+2(1 - 3r 6 ) (2/123 + 2/134) - 2/123 - 32/1232/134 - 2/134) 2/13 
+ (2 - 5r b - 22/123 - 22/134) 2/13 ~ J/13 ) 
defines the limits of the phase space. 

2. System 3-4 

For diagrams containing the denominator j/34 = 2(p 3 -p4)/m 2 ? , that is, gluons 
labelled as 3 and 4 soft and/or collinear we will work in the cm. frame of both 
gluons where 

pi = {Ex,... ,0,pi) , 

p 2 = {E 2 ,... ,p 2 sin/3,p 2 cos/3) , 

p 3 = E 3 (l, . . . , sin 8 cos 8' , cos 6) , 

Pi = Ei{l, . . . , — sin#cos#', — cos#) . 

In this case we will use the following set of variables 

2(^3 • Pi) 2(pi • p 34 ) 2(p 2 • P34) 

2/34 = 2 ' J/ 134 ~~ 2 ' V 234 ~ 2 ' 

m z m z m z ^ ^ 
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Figure 2. Threemomenta in the cm. frame of particles 3 and 4. 



with P34 = j>3 + Pi, in terms of which energies and threemomenta read 

77n 2/134 1 /~1 " 

#1 = 75-7== , Pi = 7-7-== ^134 - 4^2/34 

p 2/234 1 



E3 



VV34 



P2 



E 4 



V2/34 



2/234 - 4rf,7/34 



2 ' 2 
We obtain for the D-dimensional phase space in this system 
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16tt 2 T(l - e) 



dy 3 40(h p )h p e y 34 
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Jo N e> Jo 



with 



(B.14) 



(B.15) 



hp = 2/1342/234(1 ~ 2/134 ~ 2/234) - r fc (yi34 + 2/234) 2 

+ (-1 + 4r 5 + 2(1 - 2r 6 )(yi 34 + 2/234) - 2/134 - 32/1342/234 - 2/234) 2/34 



2(1 - 2r b - yi 34 - 2/234)2/34 - 2/34 



(B.16) 



2. SYSTEM 3-4 
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Observe h p reduce to (B.4) in the case 2/34 — > 0. Furthermore, in this limit 
P34 = P3 + P4 behaves as the momentum of a pseudo on-shell gluon because 
pIa - 0.' 



B. PHASE SPACE IN D = 4 - 2e DIMENSIONS 
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